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We report the results of a search for Z(4430)~ decay to J/ip%~ or ip(2S)n~ in — > J/ip-K~ K°' + 
and B~'° — > ip(2S)n~ K°' + decays. The data were collected with the BABAR detector at the SLAC 
PEP-II asymmetric-energy e + e~ collider operating at center of mass energy 10.58 GeV, and the 
sample corresponds to an integrated luminosity of 413 fb -1 . Each Kn~ mass distribution exhibits 
clear K*(892) and if! (1430) signals, and the efficiency-corrected spectrum is well-described by 
a superposition of the associated Breit-Wigner intensity distributions, together with an 5*-wave 
contribution obtained from the LASS 1—1/2 Ktt~ scattering amplitude measurements. Each 
K-k~ angular distribution varies significantly in structure with Kir~ mass, and is represented in 
terms of low-order Legendre polynomial moments. We find that each J/tpir~ or -4>(2S)n~ mass 
distribution is well-described by the reflection of the measured Kir" mass and angular distribution 
structures. We see no significant evidence for a Z(4430)~ signal for any of the processes investigated, 
neither in the total J/ipn" or ip(2S)n~ mass distribution, nor in the corresponding distributions 
for the regions of K-k~ mass for which observation of the Z(4430)~ signal was reported. We 
obtain branching fraction upper limits B(B~ -> ZJK°,Z' -> J/ip-w~) < 1.5 x 10~ 5 , B(B" -> 
Z~K+, Z~ -> J/V>0 < 0.4 x 10" 5 , B(B- -> Z~K°, Z~ -> V(2S>-) < 4.7 x 10" 5 , and B(B° -> 
Z~K+, Z~ -> ip(2S)ir~) < 3.1 x 10~ 5 at 95% confidence level, where the Z(4430)~ mass and width 
have been fixed to the reported central values. 
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PACS numbers: 12.39.Mk,12.40.Yx,13.25.Hw,14.40Gx 



I. INTRODUCTION 

In the original paper in which he proposed the Quark 
Model d], Gell-Mann stated that "Baryons can now be 
constructed from quarks using the combinations (qqq), 
{qqqqq), etc., while mesons are made out of (qq), (qqqq), 
etc.". He chose the lowest configurations to create the 
representations describing the known meson and baryon 
states. However, the higher configurations were not a 
priori excluded, and experimentalists and theorists have 
been seeking evidence supporting the existence of such 
states ever since. 

In the baryon sector, resonant structure in the KN 
system would be indicative of five-quark content, and 
searches for states of this type have been carried out since 
the mid-1960's, mainly through partial-wave analysis of 
KN elastic and charge-exchange scattering data. In re- 
cent years, there has been a great deal of activity focused 
on the search for the conjectured 0(154O) + pentaquark 
state decaying to K°p. However, the initial low-statistics 
signals claimed in a variety of experimental contexts have 
not withstood high-statistics scrutiny, and the existence 
of this state must be considered to be in doubt at the 
present time. The subject is reviewed in Ref. @, and the 
status is updated in Ref. [|[. 

In the meson sector, attention has been focused over 
the years mainly on the ao(980) and the /o(980) scalar 
mesons as possible four-quark states. However, the dis- 
covery of the £>* (2317) and the £> s i(2460), with their 
unexpectedly low mass values, and the observation of 
many new charmonium-like states above threshold for 
decay to open charm, have led to speculation that cer- 
tain of these may be four-quark states (see e.g. Ref. [1]), 
although in no case has this been clearly established. In 
this regard, it follows that the recent paper from the 
Belle Collaboration Q which reports the observation of a 
resonance-like structure, the Z(4430)~, in the ip(2S)ir~ 
system produced in the decays B~'° — > ip(2S)ir~ K°' + J6| 
has generated a great deal of interest (see e.g. Ref. Q, 
and references therein). Such a state must have a mini- 
mum quark content (ccdu), and would represent the un- 
equivocal manifestation of a four-quark meson state. 

It is clearly important to seek confirmation of the Belle 
observation, not only in the ip(2S)n~ system, but also for 
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the J I combination, which might also show evidence 
of a Z(4430)~ signal or of a similar lower mass state Q. 
Consequently, in this paper we present a BABAR analy- 
sis of the entire Dalitz plot corresponding to the decays 
B~'° — > ip(2S)ir~ K°' + and in parallel pursue an identical 
analysis of our B~'° — > J/i(m~ K Q:+ data. Both analy- 
ses make use of the complete BABAR data sample accrued 
at the Y(4S) resonance. In this regard, we first seek a 
representation of the Kn~ mass and angular distribution 
structures, which dominate the final states under study, 
in terms of their expected low-angular-momentum inten- 
sity contributions. We then investigate the reflection of 
each Kir~ system into its associated i{jtt~ @ mass distri- 
bution in order to establish the need for any additional 
narrow signal. 

The BABAR detector and the data sample are described 
briefly in Sec. [Til and the event selection procedures are 
discussed in Sec. IIIII In Sec. [TVl the Dalitz plots and 
their uncorrected invariant mass projections are shown 
for the B meson signal regions. Since the analysis em- 
phasizes this search for narrow structure in the J/ipir~ 
and ip(2S)ir~ mass distributions, the mass resolution de- 
pendence on invariant mass for these systems is analyzed 
in Sec. [V] Similarly, it is important to understand the 
behavior of the event reconstruction efficiency over each 
final state Dalitz plot and to correct for it before assess- 
ing the significance of any observed mass structures. The 
procedure followed is described in Appendix A, and the 
results are summarized in Sec. I VII Fits to the corrected 
Ktt~ mass distributions are discussed in Sec. IVII1 and 
the K~k~ angular distribution structure as a function of 
Ktt~ mass is represented in terms of Legendre polyno- 
mial moments as described in Sec. I Villi In Sec. IIXI the 
reflections of the observed Kir~ mass and angular struc- 
tures onto the J/ipn~ and iP(2S)tt~ mass distributions 
are compared to the corresponding efficiency-corrected 
distributions, and in Sec. [X] our results are discussed in 
relation to those in the Belle publication. The BABAR 
tpTr~ mass distributions are fitted in Sec. IXII and we 
present a summary and our conclusions in Sec. IXIII Fi- 
nally, acknowledgments are expressed in Sec. IXIIII 

II. THE BABAR DETECTOR AND DATA 
SAMPLE 

The data used in this analysis were collected with the 
BABAR detector at the PEP-II asymmetric-energy e + e~ 
collider operating at a center-of-mass (cm.) energy of 
10.58 GeV. 

A detailed description of the BaBar detector can 
be found in Ref. [10] , Charged particle tracks are 
detected with a five-layer, double-sided silicon vertex 
tracker (SVT) and a 40-layer drift chamber (DCH), filled 
with a helium-isobutane gas mixture, and coaxial with 
the cryostat of a superconducting solenoidal magnet, 
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which produces a magnetic field of approximately 1.5 T. 
The charged-particle momentum resolution is given by 
(SPt/Pt) 2 = (0.0013P T ) 2 + (0.0045) 2 , where P T is the 
transverse momentum measured in GeV/c. The SVT, 
with a typical coordinate resolution of 10 /xm, measures 
the impact parameters of charged particle tracks in both 
the plane transverse to the beam direction and along the 
collision axis; it also supports stand-alone reconstruction 
of 1ow-Pt charged particle tracks. 

Charged particle types are identified from specific ion- 
ization energy loss (dE/dx) measured in the DCH and 
SVT, and from Cherenkov radiation detected in a ring- 
imaging Cherenkov device. Electrons are identified by 
means of a CsI(Tl) electromagnetic calorimeter (EMC). 

The return yoke of the superconducting coil is instru- 
mented with resistive plate chambers for the identifica- 
tion of muons and the detection of clusters produced 
by Kl and neutron interactions. For the latter part of 
the experiment these chambers were replaced by limited 
streamer tubes in the barrel region of the detector [ll[ . 

In this analysis, we use a data sample corresponding to 
an integrated luminosity of 413 fb , which is equivalent 
to the production of approximately 455 million BB pairs. 

III. EVENT SELECTION 

We reconstruct events in four decay modes ((J: 

B~ -» J/iPtt-K", (1) 

B° -» J/iPn-K+, (2) 

B~ -> iP(2S)n-K° s , (3) 

B° -> ^{2S)ir^K+ . (4) 

The event selection criteria were established by op- 
timizing signal-to-background ratio using Monte Carlo 
(MC) simulated signal events, B~ ,a — ► tpir~K 0,+ , and 
background, BB and e + e~ — > qq (q — u, d, s, c), events. 

For the data sample, a J/ip candidate is formed by 
geometrically constraining an identified e + e~ or 
pair of tracks to a common vertex point and requir- 
ing a fit probability > 0.001. For fi + [i~, the invari- 
ant mass of the pair must in addition satisfy 3.06 < 
"V+u- < 3.14 GeV/c 2 , while for e + e~ the requirement 
is 2.95 < m e + e - < 3.14 GeV/c 2 . In the latter case, the 
mass interval extends to lower values in order to allow 
for electron bremsstrahlung energy loss; if an electron- 
associated photon cluster of this type is found in the 
EMC, its four-momentum vector is included in the calcu- 
lation of m e + e - . The surviving J/tp candidates were fit- 
ted to impose a constraint to the nominal mass value [|| . 

For ip(2S) decay to or e+e~ the same selec- 

tion procedures are followed, but with invariant mass re- 
quirements 3.640 < m < 3.740 GeV/c 2 or 3.440 < 
m e+e - < 3.740 GeV/c 2 For tp(2S) decay to J/ipn+n-, 
the J/tp candidate is selected as previously described, 
and is fit again to incorporate a constraint to its nominal 
mass value 0. This J/ip and an identified Tt + n pair are 



geometrically constrained to a common vertex (fit prob- 
ability > 0.001), and required to have an invariant mass 
in the range 3.655 < mj/^ w +„- < 3.715 GeV/c 2 . In the 
same manner as for the J/ip , surviving candidates were 
then constrained to the nominal ip(2S) mass value Q. 

A Kg candidate is formed by geometrically constrain- 
ing a pair of oppositely charged tracks to a common ver- 
tex (fit probability > 0.001); the tracks are treated as 
pions, but without particle-identification requirements, 
and the invariant mass of the pair must satisfy 0.472 < 
m„.+„— < 0.522 GeV/c 2 . A charged kaon candidate 
from the B meson decay must be identified as a kaon, 
but no particle identification is required of correspond- 
ing charged pion candidates. 

The ip, K and 7r candidates forming a B meson de- 
cay candidate are geometrically constrained to a com- 
mon vertex, with fit probability > 0.001 required. For 
decay modes involving a Kg, the K s flight length with 
respect to this vertex must have > +3 standard deviation 
significance in order to reduce combinatoric background. 
The K s candidate is not mass constrained, since this was 
found to have a negligible effect on resolution. 

We further define B meson decay candidates using 
the energy difference AE — E* B — \fs/2 in the center 
of mass (cm.) frame, and the beam-energy substituted 
mass m ES = yj ((s/2 + pi ■ p B )/E l ) 2 - pg, where (£?;,$) 
is the initial state four-momentum vector in the labora- 
tory frame and y/s is the cm. energy; E* B is the B meson 
energy in the cm. and pb is its laboratory frame mo- 
mentum. 

We require that B decay signal events satisfy 5.272 < 
toes < 5.286 GeV/c 2 and \AE\ < 0.020 GeV. In order 
to correct for background events in the signal region, we 
define a AE sideband region by 0.030 < \AE\ < 0.050 
GeV; we have verified through MC studies that sideband 
events in the toes signal range correctly represent back- 
ground in the B meson signal region. We refer to the 
procedure by which we correct for background in the sig- 
nal region by subtracting the AE sideband events in the 
Toes signal range by the term "sideband subtraction" . 

In Figs. QJa)-(d) we show the m E s distributions in the 
AE signal region for the decay processes of Eqs. JT])- Q, 
where the filled histograms show the sideband distribu- 
tions. We fit each distribution with a signal Gaussian 
function with mass and width as free parameters, and 
an ARGUS background function [12] with a free expo- 
nential slope parameter. In each figure, the solid curve 
represents the total function and the dashed curve shows 
the background contribution. Clear toes signals are ob- 
served in Figs. [Ha)-(d), and in each figure the sideband 
distribution is consistent with the fitted background. 

The AE distributions for the toes signal region 
(Figs. [He)-(h)) exhibit clear signal peaks. We fit each 
distribution with a linear background function and a sig- 
nal function consisting of two Gaussian functions with a 
common center; all parameters are free in the fits. In 
each case, the filled histogram is from the toes sideband 
region defined by 5.250 < ities < 5.264 GeV/c 2 , and is 
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FIG. 1: The ities distributions, (a)-(d), and (A_B) distributions, (e)-(h), for the decay modes B~ —* J/ipir~ K s , B° — > 
J/i/jtv~ K + , B~ — > iP(2S)tt~ Kg, and B° — > ip{2S)n~ K + . The points show the data, and the solid curves represent the 
fit functions. The dashed curves indicate the background contributions, and the filled histograms show the corresponding 
distributions for the sideband regions. 



in good agreement with the fitted background. For the 
decay modes of Eqs. (JTJ)— , the fraction of events with 
more than one B meson signal candidate ranges from 
0.5 to 1.1 %. For such events, the candidate with the 
smallest value of \AE\ is selected. 

We summarize the principal selection criteria in Ta- 
ble U and in Table |TT] provide an overview of the data 
samples in the -B-meson signal region used in the analy- 
sis described in this paper. 



IV. THE DALITZ PLOTS AND INVARIANT 
MASS PROJECTIONS 

The Dalitz plots of to 2 , _ versus m\. _ are shown in 
Fig. [5] for the signal regions defined in Table U for the B 
meson decay modes specified in Eqs. ([U-S]). The corre- 
sponding m Klx - , m^- , and m^K mass projections are 



represented by the data points in Figs. [3l 0] and [5j re- 
spectively. In each figure the filled histogram is obtained 
from the relevant AE sideband region. 

In Fig. [3J the contributions due to the if*(892) domi- 
nate the mass distributions. Small, but clear, if2(1430) 
signals are evident for the J/ip decay modes, and these 
seem to be present for the ip(2S) modes also. Previous 
analyses [l3|, [l4| have shown that, for the J/ip modes, 
the region between the if* (892) and if;? (1430) signals 
(~ 1.1 — 1.3 GeV/c 2 ) contains a significant Ktt~ S-wave 
contribution. In the if* (892) region, the presence of the 
S-w&ve amplitude has been demonstrated through its in- 
terference with the if* (892) F-wave amplitude [14] . This 
interference yields a strong forward-backward asymmetry 
in the if ir~ angular distribution, as is seen in the verti- 
cal if* (892) bands of Fig. H(b) and Fig. O^d), and as is 
shown in Sec. ECO Figs.H^a) and Fig.^c). These fea- 
tures of the if ir~ mass and angular distributions will be 
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FIG. 2: The versus Tn? Kv - Dalitz plot distributions for the signal regions for the decay modes (a) B — > J/ip-x K%. (b) 

B° -> J/ipir-K + , (c) B~ -> ?/>(2S')7r~A'g, (d) B° -> ^(2S)7r~.K + . The intensity scale is logarithmic. 



analyzed in detail in Sees. IVlJ and IVIlTl below. 

The to^tt- distributions of Fig. [4] show no peaking 
structure at the mass reported for the Z(4430)~ Q (in- 
dicated by the dashed vertical line in each figure). In 
Fig. BJb) there seems to be a peak at ~ 4.65 GeV/c 2 and 
perhaps a weaker one just below 4.4 GeV/c 2 , while in 
Figs. |3Jc) an( i (d) there seems to be a peak just below 
~ 4.5 GeV/c 2 . These features are discussed in Sees. IXl 
and lXIl in conjunction with reflections resulting from the 
K n~ mass and angular structures. 

Similarly, the m^K distributions of Fig. [5] show no 
evidence of narrow structure. In fact, in overall shape 
these distributions approximate mirror images of those 
in Fig. [4] This is not unexpected if both result primarily 
from K 7r~ reflection, since then the high-mass region of 
one distribution would be correlated strongly with the 
low mass region of the other, and vice versa. Since Fig. [5] 



shows no evidence of interesting features, and since our 
emphasis in this paper is on the search for the Z(4430)~, 
we do not discuss the tpK systems any further in the 
present analysis. 

We note that in Figs.OOIl and in other invariant mass 
distributions to follow, the J/ij)-ip(2S) mass difference 
causes significant differences in the range spanned in the 
respective decay modes. This should be kept in mind 
when making J/ip-ip{2S) comparisons. 



V. THE V"r MASS RESOLUTION 

In Ref. @, the width of the £(4430)" is given as 45±f| 
MeV, where we have combined statistical and system- 
atic errors in quadrature. This value is very similar to 
that of the if* (892) Q, although with larger uncertain- 
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TABLE I: Summary of the principal criteria used to select B 
candidates. 



Selection category 


criterion 


J/ip —> e + e~ 


2.95 < m ee < 3.14 GeV/c 2 




3.06 < m MM < 3.14 (jev/c 


tp(28) -> e+e" 


3.44 < m ee < 3.74 GeV/c 2 


tp(28) -> J/4>it + tt- 
\J/W — > e e J 


3.655 < mjyv,™ < 3.715 GeV/c 2 


ip(2S) -> fi+fj," 


3.64 < m MM < 3.74 GeV/c 2 


I ( O C\ T / / ~l~ — 

(J/iP - m+m") 


3.655 < rriji^^^ < 3.715 GeV/c 




0.472 < ro™ < 0.522 GeV/c 2 


Flight length significance 


> +3(7 


ties signal region 


5.272 < m ES < 5.286 GeV/c 2 


A_E signal region 


|A£| < 0.020 GeV 
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FIG. 4: The m^, n - mass projections for the Dalitz plots of 
Fig. [2] The data points are for the m^s-AE signal regions, 
and the filled histograms are for the AE sideband regions. 
The dashed vertical lines indicate m^- = 4.433 GeV/c 2 . 
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FIG. 3: The m KlT - mass projections for the Dalitz plots of 
Fig. The data points are for the m^s-AE signal regions, 
and the filled histograms are for the AE sideband regions. 



FIG. 5: The m^K mass projections for the Dalitz plots of 
Fig. [2] The data points are for the m^s-AE signal regions, 
and the filled histograms are for the AE sideband regions. 



ties, which we have no difficulty observing, as shown in 
Fig. [3] However, mass resolution degrades with increas- 
ing Q- value, where Q- value is the difference between the 
invariant mass value in question and the corresponding 
threshold mass value. Since the Q-value for the if* (892) 
is only ~ 260 MeV/c 2 , while that for the Z(4430)~ is 
~ 600 MeV/c 2 in the i/>(2<S> mode, and would be ~ 1200 
MeV/c 2 in the J/ifjir mode, the mass resolution should be 
worse for the Z(4430) _ than for the if* (892), and hence 
should be systematically investigated. 

We do this by using MC simulated data for the B me- 
son decay modes of Eqs. (H])-((ll). For each mode, the 
reconstructed events were divided into 50 MeV/ c 2 inter- 
vals of V'tt - mass, and within each interval it was found 
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TABLE II: The data samples used in the analysis. 



Decay mode 


Signal region events 


AE sideband events 


Net analysis sample 


B~ -> J/ij)-K-K% 


4229 ± 65 


485 ± 22 


3744 ± 68 


B° -> J/tPiv-K + 


14251 ± 119 


1269 ± 36 


12982 ± 124 


B~ -> ^(25)7T _ K° s 


703 ± 26 


161 ± 13 


542 ± 29 


B° -f 4>{2S)ti-K + 


2405 ± 49 


384 ± 20 


2021 ± 53 



that the distribution of (reconstructed - generated) ipir~ 
mass could be described by a double Gaussian function 
with a common mean at zero. The local ipn~ mass resolu- 
tion is characterized by the half-width-at-half-maximum 
(HWHM) value for this line shape, and the dependence 
of this quantity on ipn~ mass is shown in Fig. [5] for the 
individual decay modes. For both J/ip modes, the res- 
olution varies from ~ 2 MeV/c 2 at threshold to - 9 
MeV/c 2 at the maximum mass value, while for both 
ip(2S)w modes, the variation is from ~ 2 — 6 MeV/c 2 . 
At the Z(4430)~ mass value, indicated by the dashed 
vertical lines in Fig. [51 the resolution is ~ 4 MeV/c 2 for 
^(2S)n, and ~ 7 MeV/c 2 for J/^tt. We note that, for 
J/ipir, the resolution at a Q- value of ~ 600 MeV/c 2 is 
essentially the same as for iP(2S)tt at the Z(4430)~. It 
follows that failure to observe the Z(4430)~ in its J/ifjir~ 
or ip{2S)ir~ decay mode in the present experiment should 
not be attributed to inadequate mass resolution. 




3.5 4 4.5 4 4.5 

m J/w (GeV/c 2 ) m ¥2S)K (GeV/c 2 ) 

FIG. 6: The ipir~ mass resolution in 50 MeV/c 2 intervals 
as a function of ipn~ mass for the decay modes (a) B~ — > 
J/ipir'Ks, (b) B° -> J/im'K + , (c) B~ -> V(2S>~ K%, and 
(d) B° — » Tp(2S)-K~ K + . The dashed vertical lines indicate 
m^- = 4.433 GeV/c 2 . 



VI. EFFICIENCY CORRECTION 

In the search for the Z(4430)~, a detailed understand- 
ing of event reconstruction efficiency over the entire final 
state Dalitz plot for each of the B meson decay processes 
of Eqs. dU)- is necessary. This is because efficiency 
variation can, in principle, lead to the creation of spuri- 
ous signals or to the distortion of real effects such that 
their significance is reduced. Even when the process of 
efficiency correction leads to no significant change in the 
interpretation of the data, it is important to demonstrate 
clearly that this is in fact the case. 

The efficiency correction procedure which we follow is 
described in detail in Appendix lAl For reasons discussed 
there, we use a "rectangular Dalitz plot" for which the 
variables are chosen to be wik-k- and cos 9k, the nor- 
malized dot-product between the Kn~ three-momentum 
vector in the parent-B rest frame and the kaon three- 
momentum vector after a Lorentz transformation from 
the B rest frame to the Ktt~ rest frame. For the Dalitz 
plots shown in Fig. [21 the j/-axis variable, fnr,_ , varies 
linearly with cos8k- 

The average efficiency Eq depends on ip decay mode, 
as shown in Fig. [7] The individual fitted curves in this 
figure are used to calculate the value of Eq for an event at 
a particular value of m K7T - which has been reconstructed 
in the relevant ip decay mode. 

As explained in Appendix [K\ Eq is then modulated 
by a linear combination of 12 Legendre polynomials in 
cos 0k , whose multiplicative coefficients E\ — E\2 are ob- 
tained from curves representing their individual rriK-K- 
dependence. In this way, an efficiency value can be cal- 
culated for each reconstructed event in our data sam- 
ple. The inverse of this efficiency then provides a weight- 
factor which is associated with this event in any distribu- 
tion to which it contributes. This enables us to correct 
any distribution under study for efficiency-loss effects. 

A specific example is provided by the Ktt~ mass dis- 
tributions of Fig. [5] A weight value for each event is cal- 
culated according to its particular ip decay mode, as de- 
scribed above. The histograms of Fig. [5] are then formed 
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by summing these weights in each mass interval of each 
plot. Sideband subtraction is accomplished by assigning 
sideband events negative weight. The contributions from 
the different ip decay modes are distinguished by shad- 
ing, and the final histograms represent the sum of these 
contributions. 

In Appendix [X] it is pointed out that the use of 
high-order Legendre polynomials is necessary because of 
significant decrease in efficiency for cos 8k ~ +1 and 
0.72 < m KlT - < 0.92 GeV/c 2 , and for cos 6^ ~ -1 and 
0.97 < m Kw - < 1.27 GeV/c 2 . The former loss is due 
to the failure to reconstruct low-momentum pions in the 
laboratory frame, and the latter is due to a similar failure 
to reconstruct low-momentum kaons. The dependence of 
efficiency on laboratory frame momentum is shown in 
Fig. [5J In Figs.[Ua)-(d), there is a significant decrease in 
efficiency for pions of momentum < 0.1 GeV/c, while in 
Figs. IHf),(h) there is similar decrease for charged kaons 
below i 

for charged pions (Figs. EIe),(g)). 

The Lorentz boost from the laboratory frame to the 
Ktt~ rest frame translates the laboratory frame losses 
into the losses localized in rrix-n- and cos 9k which our 
efficiency study reveals. The effect of these regions of 
low efficiency on the uncorrected ipn mass distributions 
is discussed in Sec. HXl 



0.25 GeV/c. For K s , the effect is similar to that 



VII. FITS TO THE Kir~ MASS DISTRIBUTIONS 

The most striking aspects of the Dalitz plots of Fig. [2] 
pertain to the Kir~ system, as shown explicitly in Fig. [3] 
In order to investigate the features of the ipn" distribu- 



tions of Fig.|4l it is necessary to understand how structure 
in the Kir~ mass and cos 9k distributions reflects into the 
system. We begin this process by making detailed 
fits to the sideband-subtracted and efficiency-corrected 
Kir~ mass distributions of Fig. [8] 

As discussed in Sec. IIVI we expect that the Kir~ sys- 
tem can be described in terms of a superposition of S-, 
P-, and P/-wave amplitudes. Since we correct for cos 9k 
dependence of the efficiency, it follows that the corrected 
Ktt~ mass distributions of Fig. [5] can be described by 
a sum of S-, P-, and P-wave intensity contributions, 
since any interference terms vanish when integrated over 
cos 9k- Consequently, we describe the tuk-k- mass pro- 
jections as follows, 



dN 
drriK-K 



Nx 

Gs 



(5) 



/ GsAi 



j GodrriKT, 



where the integrals are over the full mj(, range and the 
fractions / are such that 



fs + fp + Id = 1 



(6) 



The P- and P-wave intensities, Gp and Gd, are ex- 
pressed in terms of the squared moduli of Breit-Wigner 
(BW) amplitudes. For the P-wave, 



Gp{m K n) 



Bp(m Kv )(p ■ q) 



Dp{qR P ) 



(7) 



(m 2 p - m 2 K J 2 + m 2 p T 2 p (q, R P ) 
where 

• Bp(rnK-n) describes the P-decay vertex; 

• p is the momentum of the tp in the B rest frame; 

• q is the momentum of the K in the Kix~ rest frame; 

• Dp is the P-wave Blatt-Weisskopf barrier factor 
with radius Rp fill; 



the mass-dependent total width is 



r P = r° « 



Dp(q P Rp) ( q \ ( m P 



qj, J Dp(qRp) \qp J \mKn 



(8) 



with qp — q evaluated at rap; rap is the mass, and T° p 
the width, of the P*(892). We leave the mass and width 
free in the fits and choose Rp — 3.0 GeV -1 fij|. 
Similarly, for the P-wave, 



B D (mK*)(jP ■ q) DD q [qRD) 
{ m 2 D -m 2 K7r ) 2 + m 2 D T 2 D (q,R D ) 

• Bp)(rriKTr) describes the P-decay vertex; 



G D (raK«) - , 2 _ :2 ., , 
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FIG. 8: The if7r mass distributions, after sideband subtraction and efficiency correction, for the decay modes (a) B — > 
J/ipn-K^, (b) B° -> J/ipw-K+, (c) B~ -> ^(25)^-^, and (d) B° -> ^(2S>-is:+. The contributions from the individual 
xj) decay modes are obtained separately, as described in the text, and are accumulated as indicated by the keys in (a) and (c), 
to form the final histograms. 



• Dd is the D-wave Blatt-Weisskopf barrier factor 
with radius Rd [3; 

• the mass-dependence of the total width is approx- 
imated by the Ktt~ contribution as follows, 

r° = M*)¥r^( ± )(—)> (10) 

\q D J D D (qR D ) \q D J \m K7T J 

with qu = q evaluated at did; mn is the mass, and rj, 
the width of the K%{ 1430). We fix m D and T° D to their 
nominal values d and choose R D = 1.5 GoV" 1 0. 

The S*-wave contribution is described using the 7=1/2 
amplitude for S'-wave K~ir + elastic scattering (l8| . We 



write 

Gs(m^) = S s (m^)(p-g)|T s | 2 , (11) 

where Bs{m,KTr) describes the -B-decay vertex; T$ is the 
invariant amplitude, which is related to As, the complex 
Kn~ scattering amplitude, by 

\Ts\ = (^f-) \A S \ • (12) 

For m K7T - > 1.5 GeV/c 2 , \Ag\ is obtained by interpo- 
lation from the measured values [Tr| . For lower mass 
values, As is a pure-elastic amplitude (within error) and 
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FIG. 9: (a)-(d): The dependence of efficiency on pion momen- 
tum in the laboratory frame for the B meson decay modes of 
Eqs. QJ-ljl]). (e)-(h): The corresponding efficiency depen- 
dence on kaon laboratory frame momentum. 



is parameterized as 
1 



A, 



2i& E 



cot Sr — i 



1 



cot Sr — i 



(13) 



where the first term is non-resonant, and the second is a 
resonant term rotated by 28b in order to maintain elastic 
unitarity. In Eq. 1131 



1 1 ^ 
gcotd B = - + -rq , 
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with a = 1.94 GeV" 1 and r = 1.76 GeV" 1 Kl 
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m s (= 1.435 GeV/c 2 ) is the mass of the KZ(U30) reso- 
nance and (= 0.279 GeV) is its width [l6fl. We choose 
to fix the S-wave parameters to the indicated values. 

Although the vertex functions Bs, Bp, and Bo de- 
pend, in principle, on mx-n, we find that our best fits 
to the mass distributions are obtained when each ver- 
tex function is set to one (Bs = Bp = Bp, = 1). In 
this regard, we emphasize that our goal is not to ob- 
tain a precise amplitude decomposition of the Kir~ mass 
spectrum. This would require taking account of the an- 
gular correlations between the tp decay products and the 



Ktt~ system. For Kir~ S-, P-, and Z?-wave this is ex- 
tremely complicated (l7| . and is far beyond the scope of 
the present analysis. Our aim is to obtain an accurate 
description of the Ktt~ mass distributions in terms of 
the expected angular momentum contributions so that 
we can reliably project the observed structures onto the 
related tpT:~ mass distributions. 

The results of the fits to the ttik-k- distributions of 
Fig. [8] are shown by the curves in Fig. [TUJ Good descrip- 
tions are obtained, even though each fit has only five 
free parameters and the fit function has exactly the same 
structure in each case. 

Table IIIII summarizes the output from the fits. The 
X 2 /NDF (NDF = Number of Degrees of Freedom) val- 
ues are satisfactory, and the mass and width values for 
the if* (892)° and K* (892) ~ modes are internally consis- 
tent and agree with their nominal values [3j; the width 
values for the ip(2S) modes are slightly low, but the un- 
certainties are quite large because of the smaller data 
sample involved. 

The fractional contributions for the two J/ip modes 
agree well with each other, and there is similar agree- 
ment for the ip(2S) modes. This demonstrates that the 
charged and neutral Kix~ distributions are very similar 
in shape, and so we combine them for the remainder of 
the analysis, unless we explicitly indicate otherwise. 

The results of repeating the fits for the combined distri- 
butions are shown in Fig. llli where we use a logarithmic 
y-axis scale in order to display the individual contribu- 
tions more clearly. The solid curves describe the distri- 
butions very well. For tuk-k ~ 0.7 GeV/c 2 , the curves in 
both plots are slightly below the data, and we believe that 
this results from the parameterization of the 5*- wave am- 
plitude at low Ktt~ mass values. If we normalize the S'- 
wave amplitudes of Ref. [jj] and Ref. [18] at 1.2 GeV/c 2 
and average them from threshold to 1.2 GeV/c 2 , the dis- 
crepancy is removed. However, the fits become slightly 
worse in the region between the two K*'s. Since the lat- 
ter region is very important to the present analysis, while 
the region around 0.7 GeV/c 2 is much less so, we do not 
make use of the modified S-wave amplitude. 

Based on Table IIIII the following branching fractions 
can be calculated: 



B{B~ -> J/iPk-R ) = (1.101 ± 0.021) x 10" 3 , 



(17) 



B(B° -> J/iPn-K+) = (1.079 ± 0.011) x 10~ 3 , 



(18) 



B(B~ -> ^(2S)n~K°) = (0.588 ± 0.034) x 10" 3 , 



(19) 

B(B° -> ip(2S)-!r-K + ) = (0.557 ± 0.016) x 10" 3 , 

(20) 

where we have corrected for the relevant J/ip_, ip(2S), 
K s — > 7r°7r°, and K\ branching fractions Q. The 
quoted errors result from the statistical uncertainties of 
Table Hill The published value corresponding to Eq. fT8")) 
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FIG. 10: The efficiency-corrected and sideband-subtracted Kn mass distributions of Fig. [5] for the decay modes (a) B — > 
J/ijiv-K^, (b) B° -> J/iP-r-K+, (c) B~ -> V(2S')7r-A'g, and (d) B° -> ip(2S)iv~ K + . The data are shown as open dots, and 
the curves correspond to the fits described in the text. 



TABLE III: Summary of the fit results corresponding to Fig. [TO] For each decay mode, we list the total number of Kn 
events after efficiency-correction and sideband-subtraction, \ 2 fNDF, A'* (892) mass and width, and the percentage S-, P-, and 
D-wave intensity contributions. Only statistical uncertainties are given. 



Mode 




Corrected Events 


X 2 /NDF 


m(A*(892)) 
(MeV/c 2 ) 


T(A*(892)) 
(MeV) 


S-wave 
(%) 


P-wave 
(%) 


D-wave 
(%) 


B~ -> 




20985±393 


117.6/149 


892.9±0.8 


49. Oil. 9 


17.0±1.6 


72.5±1.3 


10.5±1.0 


B° — > 


J/iPtv~K + 


57231±561 


171.4/149 


895.5±0.4 


48.9±1.0 


15.7±0.8 


73.5±0.7 


10.8±0.5 


B~ -> 




5016±292 


98.1/90 


891.6±2.1 


44.8±6.0 


23.4±4.5 


71.3±4.4 


5.3±2.7 


B° — > 




13237±377 


81.5/90 


895.8±1.0 


43.8±3.0 


25.4±2.2 


68.2±2.0 


6.4±1.2 



is (1.2 ± 0.6) x 10" 3 (3, and for the mode of Eq. j20]) 
an upper limit of 1 x 1CP 3 (90% confidence level (c.l.)) 
is quoted 0. No other information on these branch- 
ing fraction values exists; the present measurements thus 
represent significant improvements, even in the absence 
of systematic error studies. 

We note that the charged and neutral B meson decay 
rates to J/ipiT~K agree very well, and that this is true 
also for decay to i/j(2S)tt~ K; also, the latter decays occur 
at slightly more than half the rate of the former. 



VIII. THE Att LEGENDRE POLYNOMIAL 
MOMENTS 



At the beginning of Sec. I VIII we pointed out the need 
to understand the AT7r~ mass dependence of the angular 
structure in the Ktt~ system. In order to do this, we 
choose to represent the Af-7r~ angular distribution at a 
given m^ir- in terms of a Legendre polynomial expan- 
sion, following much the same procedure as described in 
Appendix lAl for our efficiency studies. In the notation of 
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FIG. 11: The results of the fits to the Ktv mass distributions for the combined Kir charge configurations (a) B ,0 — > 
J/ipn-K°' + and (b) B~'° -> ip(2S)ir- K°'+ . The data are shown as open dots, and the individual fit contributions are as 
indicated. 



Eq. (|A5[) . we write 
dN 



d cos 6 k 



Nj2{Pi)Pi(cos6 K ) : 



(21) 



i=0 



where N is the number of events (after correction) in a 
small mass interval centered at m K7T ~, and L = 2l max , 
where t ma x is the maximum orbital angular momentum 
required to describe the Kir~ system at m KlT - . We can 
re-write Eq. (|2"Tj) as 



d cos 6 k 



and extract the coefficients from the data using 

JV 

N{P t ) ^^2Pi{cose Kj ), (23) 
i=i 

as in Appendix [XJ We refer to this coefficient as "the 
unnormalized Pi moment" in the course of our discus- 
sion. In order to incorporate the efficiency weighting and 
sideband subtraction procedures, we extend Eq. (|23p as 
follows: 



N(Pi] 



N /1\ 

J2[-)Pi(oo86 Kj ) 



N S E 



k=l 



(24) 



where N$b is the number of sideband events falling in 
this tuk-k- interval, and the efficiency values, e, are ob- 
tained as described in Appendix [XJ For convenience, 
we introduce the notation (Pf) = N(Pi), where the su- 
perscript U indicates that we refer to the unnormalized 
moment, and re- write Eq. (|22[) as 



dN _ N 
d cos 6 k 2 



(25) 



An overview of the (jriK-n- , cos 6k) structure is pro- 
vided by the rectangular Dalitz plots (see Appendix [XJ 
of Fig.El Data for the B decay modes of Eqs. (p} and (J2]) 
have been combined in Fig. [TWa). and data for those of 
Eqs. ([3]) and ([4} have been combined in Fig. [T2T bV In 
Fig- f!2Ta) , the intensity for cos 6k < is stronger than 
for cos6 K > in the region rtiK^- < 0.85 GeV/c 2 . A 
similar asymmetry is present in the AT* (892) band, where 
in addition there is a clear decrease in intensity around 
cos 6k = 0. Above ~ 1.2 GeV/c 2 the backward region is 
enhanced, especially in the increased intensity region at 
theif^ 1430 )- For m Kir- > 1-5 GeV/c 2 the overall inten- 
sity is significantly decreased, but the backward region of 
cos 6k continues to be favored. Despite the smaller data 
sample of Fig. [T2T b). the backward region seems again 
to be favored in the if* (892) and If J (1430) regions, but 
little more can be said. 

For both Fig.[T2Ta) and[T2Tb) it is necessary to examine 
the rriK-K- dependence of the unnormalized moments in 
order to quantify these qualitative features. 
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FIG. 12: The cosOk versus m KlT - rectangular Dalitz plots for the combined decay modes (a) B~'° — > J/ipn~ K°' + , (b) 
_B~'° — > ?/;(2S)7r _ i : sr ' + . The plots are obtained after efficiency weighting, but without A_E sideband subtraction. The intensity 
scale is logarithmic, and is the same for both plots. 



In order to facilitate discussion of the rax-*- depen- 
dence of the unnormalized moments, we first express 
them in terms of S-, P- and D-wave Kir~ amplitudes. 
These expressions have been obtained from the B — > 
J/i^ttK analysis of Ref. [13], after integration over the 
J/ip decay angles; they apply equally to the B meson 
decay to i/j(2S)TvK. 

For an interval of containing N events: 



(if) 



3\l — P D cos(d Po - 6 Do ) (29) 
35 



3 V 3^^ i + 1 - D+lC0S ^ p + 1 ~ 5d +J 
P-iD-i cos(<5p_ 1 - S D _ 1 )) , 



N = Sl+PZ+Dl+Pl 1 +Pi 1 +D\ l +Di 1 (26) 2J2 

{P^) = -^Dl- 2 -^(D\ 1 + D\) . (30) 

jj The Si, Pi, and Di are amplitude magnitudes and i 

(P 1 ) = S PoCos(8 So — Sp ) (27) denotes the relevant helicity; the corresponding phase 

l~2 angles are denoted by 5 with the appropriate subscript. 

+ V 5 P ° D ° cos (° Po ~ ° Do > The helicity-zero terms, denoted in bold-face font, pro- 

/g- vide the corresponding description of K~Tt + elastic scat- 

+ \ -[P +1 D + i cos(<5p +1 — $d +1 ) tering. Equations (j26|) - (f30"|) define the five measurable 

' ^ quantities accessible to the present analysis if we restrict 

+ P-iD-i cos(Sp_ 1 - S D _ 1 )} , ourselves to S-, P-, and D-wave Kir~ amplitudes. How- 

ever, the equations involve seven amplitude magnitudes 
and six relative phase values, and so they cannot be 
solved in each m^ T - interval. For this reason, we can 
(p^ 7 ) = ,lf-p 2 _|_ v iU £} 2 ^28) only measure the mass dependence of the moments of 

" the Kir~ system, and then use the results to understand 

how the rriK-K- and cos 9k structure reflects into the ob- 
served -07T - mass distributions, as will be discussed in 
Sec. lIXI In the following, all Kir~ moments are sideband- 
subtracted and efficiency-corrected. 
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The dependence of (Pf), and of (Pg), on ffl Kr - 
is shown in Fig. [T3] For each moment, the behavior 
for tjj(2S) (Fig. IT37c).(d)) is very similar to that for 
J/tp (Fig. H3I a )j(b)), but the statistical uncertainties are 
significantly larger because the net analysis sample is 
smaller by a factor of approximately six (Table [TTJ) . 
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FIG. 13: The m fc - dependence of (a) (Pf) and (b) (Pf) 
for B° -> J/ipir-K + '°; the m x?r - dependence of (c) (Pf) 
and (d) (pf) for B°'~ -> ip(2S)iv- K+'°. 



The distributions of Fig. [TU] and Figs. [T3TT51 can be 
compared to those observed in Ref. [lq ]. Fig. 6, for the 
K~tt + system produced in the reaction K~p — > if ~7r + n 
at 11 GeV/c if~ beam momentum; the latter are rep- 
resentative of the moments structure of K~ir + elastic 
scattering. A striking overall feature is the strong sup- 
pression of the mass structure for m K7T - > 1 GeV/c 2 rel- 
ative to the if* (892) which is observed for the B meson 
decay processes of Eqs. ©-(Ill) in comparison to K~w + 
clastic scattering. 

With this in mind, the (Pf) distribution of Fig. [T3T a) 
bears a remarkable similarity to that in Ref. [l6[ , except 
that the sign is reversed. This is entirely consistent with 
the analysis of Ref. fTiJ , performed with a much smaller 
BABAR data sample (^81 fb _1 ), which showed that the 
So — Po relative phase of Eq. (f2"T)) differed by ir from that 
obtained for K~ir + elastic scattering. In the if* (892) 
region, the P-wave terms in Eq. (|2~T)) should be negli- 
gible, so that the relative phase offset should yield the 
observed sign reversal w.r.t. K~n + scattering. This sign 
reversal continues through the if! (1430) region, which 
suggests that So — Po interference remains the dominant 
contribution to (Pf), especially since {P3) is systemat- 
ically positive in this region (Figs. HHa),(c)) (Note that 
the first term in Eq. (f2"9"]l differs by only 1.8% from the 
second term in Eq. fZT} ). 



The behavior of (P2) in the if * (892) region is very 
similar to that observed in Ref. [lq . and, ignoring D- 
wave contributions, agrees in magnitude and sign with 
a calculation usin g th e values of Po, P+i, and P_i from 
Ref. [14j. In Ref. [lg|, (Pj 7 ) is positive and much larger 
at the if* (1430) than at the if * (892). Clearly, this is not 
the case in Figs.E^b) andE^d), where (Pf) is small and 
negative at the if | (1430). From Eq. (|2"g]l. this could oc- 
cur if Sq is also shifted in phase by ir relative to Po, while 
the P 2 and (D 2 +l + D 2 ^) contributions to Eq. essen- 
tially cancel. The latter is suggested by the observation 
of only a small (Pf) signal at the if |( 1430) in Fig-Efb), 
and the absence of signal in this region of Fig. [T4jd). In 
Eq. (0ni), Po is favored 3:2 over (D 2 +l + P 2 x ), whereas in 
Eq. (f2"B")) the ratio is 4:5, hence the conjecture that these 
contributions may in effect be canceling in Figs. I13f b) 
andCCJd). 




100 



-100 




100 



B"'°-> l|/(2S)Ji K + 



(d)<p>: 



1.2 1.4 1.6 1.8 2 2.2 1.2 

m Kir (GeV/c 2 ) 



1.3 1.4 1.5 1.1 
m Kir (GeV/c-) 



FIG. 14: The m Kv - dependence (for m Kw - > 1.2 GeV/c 2 ) 
of (a) (PF) and (b) (Pf) for B°'~ -*■ J/ipw~ K + '°; the 
m K7r - dependence of (c) (Pf) and (d) (Pf } for B°~ -> 



The (P3 ) and (Pj 7 ) moments require that P-wave 
contributions be present. At the statistical level of the 
present analysis, we do not expect to observe such con- 
tributions for m K7T - < 1-2 GeV/c 2 , and so we begin 
the distributions of Fig. [TJ]at this value. In Figs. HUa) 
and [T4Tc). (P3 7 ) is systematically positive for all m^.- 
mass values above 1.3 GeV/c 2 . This results from the en- 
hancement observed in Fig. [12] for cos 8k < 0, and is sim- 
ilar to the behavior in Ref. [16j . but at a much- reduced 
intensity level, as mentioned previously. The (Pj 7 ) mo- 
ments of Fig. [T4T b) and [T47 d) have been discussed al- 
ready. The dip to negative values at m^ix- ~ 1-3 GeV/c 2 
in Fig.[T4Tb) is interesting, since it may indicate that the 
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relative strength of the D 2 , and (D 2 _ 1 + D 2 ^) contribu- 
tions to Eq. (|30[) varies with m Klx - . The detailed am- 
plitude analyses of the if* (892) region [3, El do not 
consider such a possibility for the Pq, P+i, and P_i am- 
plitudes, but a mass-independent approach to these anal- 
yses would require a much larger data sample. 

Finally, the extended m K7T - range available for the 
combined B~ ,Q — » J/ifnr~ K°' + data samples allows us to 
search for evidence of F-wave amplitude contributions 
associated with the if|(1780), which has mass ~ 1.78 
GeV/c 2 and width - 0.20 GeV 0. In Fig. Qija) we 
show the m Kn - dependence of (P§) for m Kv - > 1.2 
GeV/c 2 , and in Fig. Qljb), the dependence of (P^) for 
m Kn- > 1-5 GeV/c 2 . Interference between D- and F- 
wave amplitudes could yield a (P§) distribution char- 
acterized by the underlying D-wave BW amplitude. For 
zero relative phase, the itik-k- dependence resulting from 
the overlap of the leading edge of the F-wave BW ampli- 
tude with the entire D-wave BW amplitude would resem- 
ble the real part of the D-wave BW. FigurerTSTa) exhibits 
just such behavior; the intensity increases from near zero 
to a maximum below the -fT^ (1430), passes through zero 
near the nominal mass value, reaches a minimum just 
below 1.5 GeV/c 2 , returns to zero near 1.6 GeV/c 2 , and 
has no clear structure thereafter. A (Pq) moment would 
involve Fq and (F+i + F-i) intensity contributions of 
opposite sign, just as for the P- and D-waves. We see 
no clear signal in Fig fTST b). although the distribution is 
systematically negative in the region 1.7— 1.9 GeV/c 2 , 
so it could be that these contributions almost cancel, as 
may be the case for the D-wave amplitudes. In the over- 
all rrixn- distributions, these contributions add, and it is 
interesting to note that in Fig. [TTT a) there is a small ex- 
cess of events above the fitted curve in the region of the 
K 3 (1780). Our fit to the mass distribution could pos- 
sibly be improved slightly in this region by including a 
^3 (1780) contribution, but we do not do this at present. 
It seems reasonable to interpret the vuktx- dependence of 
(P^) observed in Fig. [T5Ta) as indicating the presence of 
a ^(1780) amplitude, in which case this would be the 
first evidence for B meson decay to a final state including 
a spin three resonance. 

In summary, the angular structures observed for the 
Ktt~ systems produced in D~'° — * J/ipir~ K°' + and 
B~ : ° — > iP(2S)tt^ K°- + show interesting features, which 
can be well understood on the basis of the expected Ktt~ 
amplitude contributions. Moreover, the main features 
agree well between the J/ip and ip(2S) modes, taking 
account of the statistical limitations of the latter data 
sample. 



IX. REFLECTION OF Kn~ STRUCTURE INTO 
THE ipn- MASS DISTRIBUTIONS 

We now investigate the extent to which reflection of the 
Kir~ mass and angular structures described in Sees. IVIII 
and IVIIII are able to reproduce the efficiency-corrected 
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FIG. 15: (a) The m K ^ dependence of {P^} for m K7r - > 1.2 
GeV/c 2 for B 0, ~ -» J/tpn' K + '°; (b) the m Klt - dependence 
of (Pg) for m Kw - > 1.5 GeV/c 2 for B°-~ -> J/tPiv-K + ' . 



and sideband-subtracted ?/>7r _ mass distributions. 

We do this by using a MC generator which initially 
creates large samples of unit weight B~ Q — ► ipTT~K°' + 
events with the correct B meson production angular dis- 
tribution in the overall cm. frame, and distributed in 
Ktt~ mass according to the fit functions obtained as de- 
scribed in Sec. lVIll The distribution in cos 9k is uniform 
at each value of m K7T - , and so is described by 



dN _ N 
d cos 9 k 2 



(31) 



which is just Eq. ([25| with no angular structure. 

Since we are now dealing with efficiency-corrected dis- 
tributions, and since mass resolution cannot generate lo- 
cal structure (Fig. [6]), we need not subject the generated 
events to detector response simulation and subsequent 
event reconstruction. Consequently we can create very 
large MC samples. 

The Ktt~ mass distributions generated in this way 
for B-'° -► J/ipn-K°< + and B~'° -► V(2S>~ K°'+ are 
shown in Figs.fTBFa) and fTHT b). respectively. Each distri- 
bution contains ten million events, and has been normal- 
ized to the total number of events in the corresponding 
corrected data sample; each represents the relevant fit 
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function of Fig. [IT] very well. We can use these events to 
create mass distributions according to any desired 
selection criteria, and we provide several examples of this 
later. 
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FIG. 16: The K-k mass distributions generated according to 
the Kty~ fit functions of Fig.rfT]for (a) B''° -> J/iPtt~K°' + 
and (b) B~'° — > ip(2S)-K~ K°' + . Each distribution contains 
ten million events, normalized to the total number of events 
observed in the corrected data sample (Fig. 



However, in Sec. IVIIII we showed that there is a great 
deal of angular structure in the Kir~ system, as repre- 
sented by the m K7r - dependence of the (P/ 7 ) moments, 
and this affects the reflection from the Kit~ system into 
the tpTT~ mass distribution. 

We take this Kit~ angular structure into account by 
returning to Eq. (|25[1 . removing a factor of N/2 on the 
right side, and so obtaining 




i=l / 



i.e. 



dN 



d cos 6 k 




Pi (cos 9k) 



where 



(Pi 



2 



(32) 



(33) 



(34) 



is defined to be the "the normalized Pi moment" . We can 
thus incorporate the measured Ktt~ angular structure 
into our generator by giving weight Wj to the j event 
generated, where 



Wj 



i + J2{pf)Pi(co S e K .] 



(35) 



The (Pf) are evaluated for the m K7l - value of the j th 
event by linear interpolation of the values obtained by 
normalizing the (P/ 7 ) of Figs. 1X311X51 according to Eq. ([33)). 

The results are shown in Figs. imfTTJl where we have 
used modified mass intervals in order to reduce statisti- 
cal fluctuation. We interpolate linearly using the lines 



connecting the measured values. In order to take into 
account the statistical uncertainties, we also interpolate 
using lines connecting the +1<t error values, and lines 
connecting the — la error values, as shown by the shaded 
regions in each plot. 
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FIG. 17: The normalized moments corresponding to Fig. [13] 
obtained by using Eq. (|34[1 . Mass intervals have been com- 
bined in order to reduce statistical fluctuations. The lines 
indicate the linear interpolations used in weighting the MC 
events. The shaded regions indicate the ±1<t variations used 
to account of statistical uncertainties. 



The -07T mass distributions for the entire Kit mass 
range for the decay modes B~'° — ► J/ipir~K°' + and 
B~'° -> i;(2S)7r~K ' + are shown in Fig. H^a) and 
Fig. l20y b). respectively. The points represent the data 
after correcting for efficiency and subtracting the events 
in the AE sideband. The dashed curves show the reflec- 
tion from K ir~ assuming a flat cos Ok distribution. The 
solid curves are obtained by weighting each event accord- 
ing to Eq. (|35| . The shaded bands associated with the 
solid curves indicate the effect of interpolation using ±lu 
normalized moment values, as described above. 

We emphasize that the absolute normalization of the 
curves shown in Figs. [20Fa) and [20Fb) is established by 
the scale factor used to normalize each of our ten-million- 
event MC samples to the corresponding corrected number 
of events, as shown in Figs.fToTa'l andfTBTb'). respectively. 

Comparison of the dashed and solid curves of Fig. [20l 
shows that it is important to modulate the cos 0k distri- 
butions using the normalized moment weights of Eq. (|35| . 
Since the individual Pj (cos 0k) functions integrate to zero 
over cosOk, the incorporation of the uij weights does not 
affect the distributions of Fig. [TH This also means that 
the associated dashed and solid curves of Fig. [20] inte- 
grate to the same total number of events. 
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FIG. 18: The normalized moments corresponding to Fig. [14] 
obtained by using Eq. (|34[l . Mass intervals have been com- 
bined in order to reduce statistical fluctuations. The lines 
indicate the linear interpolations used in weighting the MC 
events. The shaded regions indicate the ±l<r variations used 
to account of statistical uncertainties. 
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FIG. 19: The normalized moments corresponding to Fig. [15] 
obtained by using Eq. (|34[1 . Mass intervals have been com- 
bined in order to reduce statistical fluctuations. The lines 
indicate the linear interpolations used in weighting the MC 
events. The shaded regions indicate the ±l<r variations used 
to account of statistical uncertainties. 



The reasons for the enhancement of the solid curves at 
high ipn" mass values, and their suppression at lower 
mass values are made clear by the rectangular Dalitz 
plots of Figs. l2~TT a) and l2~TT b). We plot cos 9^, against 
m Tp-rr- j where cos 6$ is the normalized dot-product of the 
ipTr~ three-momentum vector in the parent B meson rest 
frame and the ip three-momentum vector in the rest 
frame. A guide to the structures observed in these plots is 
provided by Figs. [2"TT c) and [2TT d). where we indicate the 
locus of the A*(892) band (0.795 - 0.995 GeV/c 2 ) and 
that of the #5(1430) band (1.330 - 1.530 GeV/c 2 ), as 
chosen in Ref. [5[; in addition, we label regions (A)-(E), 
defined in Sec. [X] The dashed vertical lines indicate the 
Z(4430)~ region from Ref. @. The high-m^^- region of 
each of these bands corresponds to cos 9k < 0, and in 
Fig. l2"TTa) this region is clearly populated preferentially 
for both K* bands, and even for the Kn~ mass range 
in between. This corresponds to the backward-forward 
asymmetry observed in Fig. 1121 and to the negative val- 
ues observed for (Pf 7 ) in Fig. H5T a). and to the positive 
values of (P^) in Fig.fHTa). These high-m^- enhance- 
ments are compensated by the low-m^- suppression of 
the solid curve relative to the dashed curve in Fig. l2"0T a) , 
since the integral along any m^- locus is independent of 
the cos 9k distribution. Similar behavior is observed for 
Fig. I2"07 b). but at a reduced statistical level. Backward- 
forward asymmetry is observed in Fig.[T2Tb), where (P| ) 
is primarily negative in Fig. ll3f c) and (Pj 7 ) is positive in 
Fig. flUf c). so that the net effect on the ipTT~ mass distri- 
bution of Fig. l2"0Tb) is much the same as in Fig. l2"07 aV 

In Figs. l2"0Tc) and l2"0Td). we show the residuals (data 
- solid curves) for the distributions in Figs. I2"07 a) and 
l2"07 b). respectively. The dashed vertical lines indicate 
m^ w - = 4.433 GeV/c 2 [|. 

There is an excess of events in Fig. [20T c) for mj/^ w - 
~ 4.61 GeV/c 2 , as well as in Fig.^b) and Fig.^c) be- 
low. The effect is associated with the if* (892) region of 
A'7r~ mass (Fig. l2~4Tb)), for which the mji^- distribu- 
tion decreases steeply in the high-mass region. The mass 
resolution there is ~ 9 MeV/c 2 (Fig. E)Ja),(c)) and we 
have not incorporated this into our calculations. For this 
reason, the calculated curve falls systematically below 
the data, hence yielding a spurious peak in the residual 
distribution. For the corresponding tn^^S)-*- distribu- 
tions, the mass resolution is ~ 4 MeV/c 2 (Fig. [6](b),(d)) 
at ~ 4.6 GeV/c 2 , the statistical fluctuations are larger, 
and no similar effect is observed (Figs. l2"GTd) . l2"5Tg) . and 
l2"6T d)). Apart from this, the distribution of the residuals 
shows no evidence of statistically significant departure 
from zero at any J/^>tt~ mass value. 

In Fig. [20Tb'). and correspondingly in Fig. [20Td). the 
small excess of events at ~ 4.48 GeV/c 2 provides the only 
indication of a narrow signal. As shown in Sec. IXI Bl this 
yields a 2.7a enhancement with mass ~ 4.476 GeV/c 2 
and width consistent with that reported in Ref. Q. The 
dot-dashed curve in Fig. I2"07 b) was obtained from the 
dashed curve by modulating the Kir~ angular distribu- 
tion using instead the normalized Kn~ moments from 
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FIG. 20: The ^vr" mass distributions for the combined decay modes (a) B~-° -> J/-07T" and (b) -> V(2S)tt"A: ' + . 

The points show the data after efficiency correction and AE sideband subtraction. The dashed curves show the Kir~ reflection 
for a flat cos 9k distribution, while the solid curves show the result of cos 9k weighting. The shaded bands represent the effect 
of statistical uncertainty on the normalized moments. In (b), the dot-dashed curve indicates the effect of weighting with the 
normalized J/ipw~K moments. The dashed vertical lines indicate the value of m^ n - = 4.433 GeV/c 2 . In (c) and (d), we show 
the residuals (data-solid curve) for (a) and (b), respectively. 



the — > J/ipir~K°' + data, which show no evidence 
of a Z(4430)~ signal. This curve and the solid curve differ 
only slightly in the range ~ 4.2 GeV/c 2 to ~ 4.55 GeV/c 2 , 
so that the Ktt~ background function at ~ 4.48 GeV/c 2 
is not very sensitive to the modulation procedure, nor to 
the presence of a small, narrow WW2SW- enhancement 
(see Sec. lXIBl for a quantitative discussion). 

We conclude that the to^msw- distribution of 
Fig. l2"0Tb). and the residual distribution of Fig. l2"D7 d). 
do not provide confirmation of the Z(4430)~ signal re- 
ported in Ref. [5j. 



X. COMPARISON TO THE BELLE RESULTS 

We now compare our results to those obtained by Belle 
for B -> ^{2S)tt-K @. 



A. The ip"K mass resolution 

In Sec. El we showed (Fig. [5]) our mass resolution de- 
pendence on Q-value, and obtained HWHM ~ 4 MeV/ c 2 
for the V(2S , )tt- system at the Z(4430)~. In Ref. [|, 
it is stated only that the mass resolution is 2.5 MeV/c 2 . 
Since the width of the Z(4430)~ is - 45 MeV 0], mass 



resolution should not be an issue for the comparison of 
similar data samples (see Sec. IX E| ). 

B. Efficiency 

We have made a detailed study of efficiency over each 
Dalitz plot for each J/ip and ip(2S) decay mode sepa- 
rately (Sec. IVip . and have identified efficiency losses as- 
sociated with low-momentum pions and kaons in the lab- 
oratory frame (Fig. We illustrate the effect of such 
losses on the m^ v - distributions using our ten-million- 
event B~'° — > ipTT~K°' + samples weighted to take ac- 
count of the Kir~ angular structure (Sec.HX|). In Fig. |2"21 
we show the ipir~~ distributions obtained as for Fig. |2"01 
(solid curves). We then require that the momentum of 
the 7T be less than 100 MeV/c in the laboratory frame 
(Fig. [5]) and obtain the shaded distributions in the m^- 
threshold regions. Similarly, the requirement that the 
kaon momentum be less than 250 MeV/c in the labora- 
tory frame (Fig. \§§ yields the cross-hatched regions near 
maximum ipTt~ mass [20(. It follows that the regions of 
lower efficiency discussed in Appendix [X] should have no 
significant effect on the region of the Z(4430)~. 

As a direct check of the effect of our efficiency- 
correction procedure, we show our m^- distribu- 
tions before and after correction in Figs. l2"3Ta),(b) 



23 




FIG. 21: The cosfl^ versus m^- rectangular Dalitz plots for (a) B~'° -> J/^tt" K°' + , and (b) -> ip(2S)-n-- K°' + ; (c) 

and (d), the corresponding plots indicating the loci of the _K"*(892) and K^iUSO) resonance bands defined in the text; regions 
A-E, defined by Eqs. (|36[) - (|40[ ). are indicated. The dashed vertical lines show the mass range 4.400 < TOj, t - < 4.460 GeV/c 2 . 
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FIG. 22: (a) The curve of Fig. \M.&), and (b) that of 
Fig. l2"0T b). obtained by cos 6k weighting; the shaded regions 
near threshold correspond to P(ir~) < 0.1 GeV/c in the lab- 
oratory frame, while the cross-hatched regions near maxi- 
mum i/j7r _ mass represent P(K) < 0.25 GeV/c in the labora- 
tory frame; the dashed vertical lines indicate m^- — 4.433 
GeV/c 2 . 



and i(d),(e) for B~>° -> J/ipn-K ^ and B->° -> 
^(25)^-^°'+, respectively. Then in Fig. [23£c) and^f) 
we show the ratio of the uncorrected and corrected distri- 
butions as a measure of average efficiency. As expected, 
the average value decreases rapidly near threshold, and 
near the maximum value for both distributions. Away 
from these regions, the efficiency increases slowly with 
increasing mass. We conclude that our event reconstruc- 
tion efficiency should have no effect on any Z(4430)~ 
signal in our data. 



In Fig.l2"37f) we show the result of a linear fit, excluding 
the regions m.^^s)-*- < 3-9 GeV/c 2 and m^^s)-*- > 4.71 
GeV/c 2 , which are seriously affected by the loss of low 
momentum pions and kaons, respectively. The fitted effi- 
ciency value increases from 13.7 % to 15.0 % over the 
fitted region. The low-efficiency regions are excluded 



when we compare our uncorrected m 



distribution 



to that from Belle (Sec. IX Ep . This is due to the fact 
that for both experiments the reconstruction efficiency 
for very low momentum charged-particle tracks in the 
laboratory frame decreases rapidly to zero (c/. Fig. [5]). 
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There is no detailed discussion of efficiency in Ref. Q. 



D. Regions of Kn mass 
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FIG. 23: The cumulative i/jtt~ mass distributions for the K s 
and if + decay modes (a) and (d) before, and (b) and (e) after, 
efficiency correction. The ratio of uncorrected to corrected 
data is shown in (c) and (f ) . The fitted line in (f ) is described 
in the text. 



C. The Kn mass and cos 9k structures 

In Ref. H, the if* (892) and 2^(1430) regions 
of the Dalitz plot are removed, and the remain- 
ing non-.Z(4430)~ contribution to the to^,( 2 5) 7T - mass 
distribution is described by a second-order polyno- 
mial in m^s)^- multiplied by the momentum of the 
ijj(2S) in the iP(2S)tt~ rest frame and by the factor 
Jm-max — Ti^(2S)w- i f° r two-out-of-three phase space, 
this latter factor should be the momentum of the recoil 
K in the B rest frame. Also, the AE sideband contri- 
butions are not subtracted prior to their fit. There is no 
discussion of a possible description of the background in 
terms of Kn~ mass and angular structures. 

In our analysis, we have considered the effect of the 
Ktt~ mass and cosf?^ structures (Sec. lVIll and Sec. I Villi) , 
and have shown how the observed features affect the 
i/jtt~ mass distributions in Sec. IIXI The resulting dashed 
curves of Fig. |2"21 cannot be described in terms of second- 
order polynomials. However, each corresponds to a pro- 
jection of the entire Dalitz plot. In order to make a direct 
comparison to the Belle data, we investigate the relevant 
regions of Kir~ mass in the following section. 



The ip(2S)n mass distribution shown in Fig. 2 of Ref. Q 
has a "K* veto" applied. This means that events within 
100 MeV/c 2 of the if* (892) or the if 2 *(1430) have been 
removed, and hence that the Kix~ mass range has, in 
effect, been divided into five regions, as follows: 



region A 
region B 
region C 
region D 
region E 



m K7r - < 0.795 GeV/c 2 , 
0.795 < m Kw - < 0.995 GeV/c 2 . 
0.995 < rriKK- < 1.332 GeV/c 2 . 
1.332 < m K T<- < 1.532 GeV/c 2 . 
m K7r - > 1.532 GeV/c 2 . 



(36) 
(37) 
(38) 
(39) 
(40) 



These regions are labeled in Fig. [2lT c) and Fig. I2l7 d). 
The ip(2S) mass distribution of Ref. [5j thus contains 
events (with sideband contribution) from regions A, C, 
and E. 

In Fig. rjH we show the corrected m^- distributions 
for regions A-E of Ktt~ mass. The solid curves and 
shaded bands correspond to those in Fig. [20j with the 
same overall normalization constants as obtained for 
Fig. 1161 i.e., there is no renormalization in the separate 
Ktt~ mass regions. In Figs. [24jf)-(j), the dot-dashed 
curves were obtained using the normalized moments from 
B~'° — > J/ipn~K°' + in conjunction with Eq. (f3"5"|) . For 
regions A and B, there is almost no difference between 
the solid and dot-dashed curves, while in the other re- 
gions the differences are less than, or of the order of, the 
statistical fluctuations in the associated data. The resid- 
uals (obtained by subtracting the solid curves from the 
data) corresponding to Fig. [24] are shown in Fig. [25] and 
show no evidence of structure. In Figs. [26] and Fig. [27] 
we make similar comparisons for the combined data in 
the K* regions (B and D), and for the if*-veto region 
(A, C, and E). Again the residuals reveal no significant 
structure. 



E. Direct comparison 

The ip(2S)iT~ mass distribution of Fig. l2"5Ta) is a re- 
production of that in Ref. except for the addition 
of error bars [2l[. In Fig. l2"5Tb) we show the equiva- 
lent distribution for our combined analysis samples for 
the B meson decay processes of Eqs. §5§ and ^ for the 
K*-veto region (A, C, and E). The mass intervals are 
the same as for Fig. l2"87 a). but no efficiency-correction 
has been performed. As mentioned in Sec. IX C[ when 
we make quantitative comparisons between Fig. I28f a) 
and Fig. I28f b) we exclude the low-efficiency regions near 
threshold and at high mass, and use only the region 
3.9 < to 



ip(2S)n- < 4-71 GeV/c 2 . We make a global com- 
parison of the data samples in Table HVl The BABAR 
sample contains ~ 8% more background than does the 
Belle sample. The net signal ratio is 1.18 ± 0.09 in favor 
of Belle, although the corresponding integrated luminos- 
ity ratio is 1.46. It follows that for BABAR, net signal per 
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FIG. 24: The ip-K~ mass distributions in regions A-E of K-k~ mass for the combined decay modes (a-e) B~~'° — > J /ij}ir~ K°' + , 
and (f-j) _B _ '° — > i/j(2S)it~ K°' + ; the open dots represent the data, and the solid curves and shaded bands are as in Fig. I20I 
but calculated for the relevant Ktt~ mass region using the same overall normalization constant as for Fig. I16I In (f-j), the 
dot-dashed curves are obtained using Ktv~ normalized moments for B~'° — ■> J/ipiv~ K°' + in Eq. (|35[) , instead of those from 
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FIG. 25: The residuals (data - solid curve) corresponding to Fig. I24I the dot-dashed vertical lines indicate m^- = 4.433 
GeV/c 2 . 
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FIG. 26: The ipn~ mass distributions for (a) B~'° -> J/^n~K°' + , and (b) B~ -> V'(2S)7r"i4' ' + , for m K7r - regions B and 
D combined; the open dots represent the data, the solid, dashed, and dot-dashed curves, and the shaded bands, correspond to 

= 4.433 



those of Fig. I20f a').(b'); (c) and (d), the corresponding residual distributions. The dashed vertical lines indicate 
GeV/c 2 . 




FIG. 27: The Tpiv' mass distributions for (a) B~ ° J/ipTv~ K°' + , and (b) B~'° -> i/j(2S)ir~ K° !+ , for m Kv - regions A, C, and 
E, combined; the open dots represent the data, the solid, dashed, and dot-dashed curves, and the shaded bands, correspond to 
those of Fig. [23 a ) j (b); (c) and (d), the corresponding residual distributions. The dashed vertical lines indicate m^ v - = 4.433 
GeV/c 2 . 
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unit luminosity is ~ 1.34, while for Belle it is 1.08, and 
that this significant increase in signal yield comes at the 
cost of only a modest increase in background level. For 
both experiments the background distribution increases 
toward threshold, and differs markedly in ip(2S)iT~ mass 
dependence from the signal. 
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FIG. 28: a) The ip{2S)n~ mass distribution after K* veto 
from Ref. [5|; the data points represent the signal region 
(we have assigned \/~N errors at present), and the shaded 
histogram represents the background contribution estimated 
from the AE sideband regions; b) shows the corresponding 
distribution from the BABAR analysis. The dashed vertical 
line indicates TW2SW- = 4.433 GeV/c 2 . 



We have shown our efficiency-corrected and sideband- 
subtracted mass distribution corresponding to Fig. 128( b) 
in Fig. HTJb), and should compare the latter to the 
equivalent Belle distribution. However, this is not avail- 
able, and so we make do with the distributions of 
Fig. [28] instead. In order to justify the use of \/N er- 
ror assignments, we combine adjacent mass intervals for 
"V(2S)tt- < 4.18 GeV/c 2 and m^ {2S)7r - > 4.55 GeV/c 2 
so that we obtain at least ten events (signal + sideband) 
in each mass interval. We then create the sideband- 
subtracted distributions of Fig. a) and Fig.l2"97b) from 
the data of Fig. 1281 a) and Fig. 1257 b). respectively. In 
Fig. I2"97 b). we have scaled our data by the factor 1.18 to 
compensate for the statistical difference between the ex- 
periments. In Fig. l297 c) we show the result of subtracting 
the distribution of Fig. [^b) from that of Fig.l2"9Ta). with 
errors combined in quadrature. There is no evidence of 
any statistically significant difference, in particular near 
m ii(2S)n- — 4.433 GeV/c 2 , indicated by the dashed ver- 
tical line. If the low-efficiency (Fig. l2"3T f)) cross-hatched 
regions of Fig. I2"97 c) are excluded, the \ 2 ~ value for the 
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FIG. 29: a) The distribution of Fig. l28T a) after combining 
mass intervals for m^^ 2 s)w- < 4.18 GeV/c 2 and ttim^sW- > 
4.55 GeV/c 2 , and carrying out sideband subtraction, b) The 
distribution of Fig. I28f b) after following the same procedure, 
and in addition scaling by 1.18, as described in the text, c) 
The difference, (a)-(b), where the errors have been combined 
in quadrature; x 2 /NDF = 54.7/58 (Probability=59.9%) ex- 
cluding the low-efficiency regions (cross-hatched). 



remaining region is found to be 54.7 for 59 mass inter- 
vals. There is one normalization constant (1.18), so that 
the comparison yields \ 2 /NDF — 54.7/58, with corre- 
sponding probability 59.9%. 

We conclude that the Belle and BABAR distributions 
of Fig. [29] are statistically consistent. We have shown 
in Fig. [20] and Figs. 1211271 that all of our corrected A" - 
distributions are well-described by reflections of the mass 
and angular structures of the Kn~ system. We refer to 
this as our Kir~ background, and in Sec. IXII we quan- 
tify the extent to which an additional Z(4430)~ signal 
is required to describe the corrected BABAR A - mass 
distributions. 

We have mentioned previously that it is the backward- 
forward asymmetry in the Kn~ angular distribution as 
a function of tcik-k- which yields the high mass enhance- 
ments seen in our ipit~~ mass distributions. We show this 
effect explicitly in Fig. [30] where we plot the distribution 
of cos On — — cos 9k for regions A, C, and E of Ktt~ mass. 
For cos^tt ~ 1, m^^sw- is near its maximum value, 
and is related linearly to cos 6^, and so it is not 

surprising that the sideband subtracted distributions of 
Fig.[307b) and Fig.[30Td) bear a strong shape resemblance 
to the corresponding rn^^s)^- distributions of Fig. [27T a) 
and Fig. l2"7T b). respectively. We note also that the in- 
crease in the sideband distribution of Fig. I30f c) towards 
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TABLE IV: Summary of m^^Sjir- data from regions A, C, and E, combined (Fig. I28[) for Belle and BABAR in the mass range 
3.9-4.71 GeV/c 2 . 

Category Belle [Fig. [M.&)] BABAR [Fig. [28Tbl] 

Total signal region events (N) 824 ± 29 786 ± 28 

Sideband contribution (B) 172 ± 9 234 ± 15 

B/N 21.5% 29.8% 

Net signal 652 ± 30 552 ± 32 



cos0,r ~ — 1 corresponds to the increase in Fig. 130Tb) 
toward ipn" threshold. Figure [501 illustrates how impor- 
tant it is to take into account the angular structure in 
the Ktt~ system in creating the shape of the associated 
i\)-k~ mass distribution, even after the removal of the K* 
regions. 
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FIG. 30: The cos#,r distributions for the combined Kir mass 
regions A, C, and E, (a), (b) for the decay modes B~'° — > 
J/iIjtv-K ' + , and (c), (d) for B~'° -> i}}{2S)-k- K°'+ . In (a) 
and (c), the points represent the data in the signal region, 
and the shaded histograms show the AE sideband events. In 
(b) and (d), the signal region distributions are shown after 
sideband subtraction. No efficiency correction is applied to 
the data. 



XI. FITS TO THE CORRECTED Vtt MASS 
DISTRIBUTIONS 



A. The fit results 

In Fig. [31] we show the results of the fits to the cor- 
rected BABAR data of Figs.[Ha),(b), Figs.|HKa),(b), and 
Figs. l2"TT a).(bh In each fit, the Ktt~ background distri- 
bution has been multiplied by a free normalization pa- 
rameter; an S-wave BW line shape, with free normaliza- 
tion, mass and width parameters, has been added inco- 
herently in order to quantify the search for a Z(4430)~ 
signal. In each figure, the solid curve represents the 
fit result; the parameter values for the corresponding 
Z(4430)~ signal are summarized in Table IVl Table fVl 
also contains the results obtained when the mass of the 
£(4430)- is fixed (4.433 GeV/c 2 ), and when both the 
mass and width (45 MeV) are fixed For all fits, 
X 2 /NDF is acceptable, but deteriorates slightly, or fails 
to improve, as first the mass is fixed and then the mass 
and width are fixed. We start each fit at the mass and 
width values of m = 4.433 GeV/c 2 and r = 45 MeV. 

For the fit to the J/i/'tt - mass distribution using the 
i4T*-veto sample with mass and width free, the obtained 
mass is ~ 100 MeV/c 2 larger than that of the Z(4430)~, 
the width is essentially undetermined, the signal is more 
than 2cr negative, and it remains at least la negative for 
the other fits. 

For the K* region, the fitted mass value is closer to 
that of the Z(4430)~, but otherwise the results are very 
similar to those for the K*-veto region in that negative 
signal values are obtained. 

For the total sample, the results are no better. The 
signal is negative by ~ 2.1cr, and remains negative by at 
least 2a as the constraints are applied. 

We conclude that there is no evidence for 
Z(4430)~ production via the decay sequence 
B~'° -> Z{U3Q)-K°'+, Z(4430)~-> J/^tt". 

For the fit to the ip(2S)ir~ K*-veto sample (which is 
equivalent to the Belle analysis sample, but sideband- 
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FIG. 31: The results of the fits to the corrected mass distributions, (a)-(c) for J/ipTr~ , and (d)-(f) for ip(2S)n~ . In each 
figure, the open dots represent the data, and the solid curve represents the fit function, which consists of the Kn~ background 
distribution with normalization free, and a relativistic S-wave BW line shape, with mass, width, and normalization free; the 
dotted curves indicate the Kir~ background functions; the dashed vertical lines indicate m^ n ~ — 4.433 GeV/c 2 . 



subtracted and efficiency-corrected) , we obtain mass and 
width values which are consistent with theirs, but a pos- 
itive signal which is only ~ 1.9a from zero; fixing mass 
and width increases this to only ~ 3.1tr. Since our ef- 
ficiency in the Z(4430)~ region is almost constant (cf. 
Fig. l25T f)). our corrected signal size with mass and width 
free (426 ± 229 events) corresponds to ^ 61 ± 33 ob- 
served events. This converts to a signal of ~ 70 events in 
Fig. 129T b). and we estimate a similar value for the Belle 
distribution of Fig. [29T a). The signal size reported in 
Ref. [5| is ~ 120 events, obtained on the basis of a back- 
ground description which ignores Kir~ mass and angular 
structures. It is interesting that we find a small posi- 
tive signal with mass and width consistent with the Belle 
values. However, with the present small data sample, 
it seems impossible to decide whether this is due to the 
production of a real state, or to the imprecision of the 
normalized Kn~ moments, primarily in the region be- 
tween the two K*'s. 

For the K* region, our fitted mass value is ~ 50 
McV/c 2 higher than the Belle value, and the signal de- 
viates from zero only by ~ 2.5a; the imposition of Belle 
mass and width values yields a signal which is less than 



~ 0.5cr from zero. 

We note that each of these regions corresponds to ap- 
proximately half of the cosO^^s) range (cf. Fig. I2"l7d)), 
so that for a flat Z(4430)~ angular distribution, we would 
expect naively that the signal size for each region would 
be the same. The values in Table [V] are consistent with 
this, however the central mass values differ by ~ 5 stan- 
dard deviations. Although there could be significant in- 
terference effects associated with the K* regions, it does 
not seem possible that these could cause the signal to 
be displaced by approximately one full-width. This, to- 
gether with the fact that both signals are in the 2-3 stan- 
dard deviation range, suggests that a more likely inter- 
pretation is that both are simply statistical fluctuations. 

Finally, for the complete sample, our fitted mass value 
is ~ 40 MeV/c 2 higher than the Belle value; the width 
is consistent with Belle's, but the signal size is only ~ 
2.7a from zero; this is reduced to ~ 2.1cr when the Belle 
parameters are imposed. 

We conclude that there is no convincing evidence for 
production of the Z(4430)~ state via the decay sequence 
B~'° -> Z(U30)-K^+, Z(4430)~->^(25)tt-, especially 
since the fits take no account of the Kn~ background un- 
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certainty associated with the normalized Kir moments. 



B. Can the Kn background absorb a Z(4430) 
signal? 

Since our evidence for the existence of a Z(4430)~ sig- 
nal in the B~'° — > ip(2S)Tr~ K 0,+ data sample is less than 
compelling, it is reasonable to ask whether our use of 
the normalized Kn~ moments to modulate the iP(2S)tt~ 
background shape might have removed part, or all, of a 
real Z(4430)~ signal. 

Firstly, our Ktt~ background curves in Figs. [501 123 
and [27] show no tendency to peak in a narrow region 
around the Z(4430) _ mass position, indicated by the 
dashed vertical line in each figure. Secondly, for the 
i/j(2S)tt~ distributions in these figures, the dot-dashed 
curves obtained using the normalized Ktt~ moments 
from the B~ ° — > J/ipir~K 0,+ modes, for which there is 
no Z(4430)~ signal, do not differ significantly in shape in 
the signal region from the solid curves from the B~>° — > 
iP{2S)tt~K°' + modes. Finally, in Ref. @ it is stated that 
it is not possible to create a narrow peak from the S-,P- 
, and £>-wave amplitude structure of the Kn~ system. 
We agree with this statement and provide a quantita- 
tive demonstration below. However, we first show how 
a Z(4430)~ signal would affect the cos8k versus rriK V - 
Dalitz plot. 

To this end, we have generated a MC sample of events 
corresponding to B~'° -> Z(4430)~Jf °'+, Z(4430)~ -> 
t/)(2S)w~ where the Z(4430)~ has the Belle central mass 
and width values, and decays isotropically. Figure I3"2l 
shows the Dalitz plot which results. The Z(4430)~ events 
yield a narrow locus confined almost entirely to the region 
cos9k < 0. 

For the Legendre polynomials which we have used to 
modulate the Ktt~ angular distribution, such a distribu- 
tion yields the following behavior: 

• Pi(cos 9k) is almost always negative; 

• P2(cos9k), which is negative for |cos#/f| < 0.58, 
is negative for ~ 0.75 < thk-k < 1-55 GeV/c 2 i.e. 
almost always; 

• P3(cos9k), which is positive for —0.78 < cos 9k < 
0, is positive over almost the entire n%K n - range 
above 1.2 GeV/c 2 , which is where we make use of 
P3(cos9k) in the data; 

• P^cos 9k), which is positive for | cos 8k\ < 0.33 
and | cos8k\ > 0.88, is mainly positive for m K7r - > 
1.2 GeV/c 2 . 

We estimate the effect of such a Z(4430)~ signal on our 
Ktt~ background by adding a MC-generated Z(4430)~ 
signal incoherently to the corrected ip(2S)ir~ data of 
Fig. [217Tb'). and subjecting this new sample to the pro- 
cedure for creating our Ktt~ background contribution to 




FIG. 32: The cos 8k versus m Kn - rectangular Dalitz plot for 
the -» Z(4430)"A"°' + , Z(4430)~-> i)(2S)w~ MC events 
generated using mz = 4.433 GeV/c 2 , Y z = 0.045 GeV, and 
with isotropic Z(4430)~ decay. 



the ijj(2S)ix mass distribution. As for Fig. [321 we use 
m z = 4433 MeV/c 2 and T z = 45 MeV in the simulation, 
and generate a flat cos 9^ angular distribution. On the 
basis of the Belle result, we would expect an observed sig- 
nal of ~ 200 events for the full cos 9$ angular range (c/. 
Fig-EI^d)), and this would yield an efficiency-corrected 
signal of ~ 1500 events (cf. Fig. [23Tf)l. Consequently, 
we generate 1500 events, of which 1493 events are within 
the (mESj AS) signal region. The m^^s)*- distribu- 
tion for these events is shown as the shaded histogram 
of Fig. [33T ah the points with error bars show the re- 
sult of combining these MC events with our corrected 
ip{2S)ir- data distribution (Fig. [20Tb)). The m fo - dis- 
tribution for this combined sample is shown in Fig. 1331 b). 
where the solid curve shows the result of a fit to the com- 
bined data, and the shaded histogram indicates the re- 
flection of the Z(4430)~ MC events. We use this new 
fit curve to generate 10 million events corresponding to 
B~ ° — > ip(2S)ir~ K°' + as before, and normalize this sam- 
ple to the combined data and MC-generated Z(4430)~ 
sample. We add the unnormalized Kir~ moments for the 
simulated Z(4430)~ events to the corrected Kw~ data 
moments, and use the distribution of Fig. l337 b) to create 
new normalized Kir~ moments as in Sec. 11X1 Finally, we 
follow the weighting procedure described in Sec. [IX] to 
create the new Ktt~ background description shown by 
the dashed curve of Fig. [34] in comparison to the com- 
bined to^(2S)7t- distribution of Fig. [337 a). 

Clearly, this dashed curve does not describe the 
Z(4430)~ signal region. However, the shape of the curve 
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TABLE V: Results of the fits to the corrected BABAR ifjir mass distributions; all fits use the relevant BABAR Kn background 
shape. 



Type of fit 


total 


K* region 


K* veto 


J/i/jtt~ sample; 

Z(4430)~ mass and width free 


m = 4455±8 MeV/c 2 
T = 42 ± 27 MeV 
N z = -901 ± 420 
X 2 /NDF = 131/154 


m = 4454 ± 4 MeV/c 2 
r = 17 ± 12 MeV 
N z = -514 ± 223 
X 2 /NDF = 135/154 


m = 4545 ± 30 MeV/c 2 
r = 100 ± 96 MeV 
N z = -411 ± 181 
X 2 /NDF = 159/154 


J/ipir~ sample; 

z?(4430) - fixed mass and free width 


r = 82 ± 77 MeV 
N z = -1098 ±490 
X 2 /NDF = 136/155 


T = 79 ± 39 MeV 
N z = -810 ± 393 
X 2 /NDF = 143/155 


T = 10 ± 12 MeV 
N z = -86 ± 64 
X 2 /NDF = 162/155 


J/ipir~ sample; 

Z(4430)~ fixed mass and width 


N z = -704 ± 249 
X 2 /NDF = 137/156 


N z = -540 ± 225 
X 2 /NDF = 144/156 


N z = -147 ± 105 
X 2 /NDF = 164/156 


■4>(2S)n~ sample; 

Z(4430) _ mass and width free 


m = 4476±8 MeV/c 2 
T = 32 ± 16 MeV 
N z = 703 ± 260 
X 2 /NDF = 93/96 


m = 4483±3 MeV/c 2 
r = 17 ± 12 MeV 
N z = 447 ± 177 
X 2 /NDF = 91/96 


m = 4439 ± 8 MeV/c 2 
T = 41 ± 33 MeV 
N z = 426 ± 229 
X 2 /NDF = 106/96 


■4>(2S)tt~ sample; 

Z(4430) ~ fixed mass and free width 


T = 97 ± 77 MeV 
N z = 710 ± 440 
X 2 /NDF = 101/97 


r = 100 ± 82 MeV 

N z = 246 ± 247 
X 2 /NDF = 102/97 


r = 36 ± 26 MeV 
N z = 414 ± 194 
X 2 /NDF = 107/97 


■4>(2S)n~ sample; 

Z(4430)~ fixed mass and width 


N z = 440 ± 212 
X 2 /NDF = 101/98 


N z = 89 ± 162 
X 2 /NDF = 101/98 


N z = 431 ± 137 
X 2 /NDF = 107/98 



has been changed compared to that of Fig. [20f b) be- 
cause of the effect of the Z(4430)~ signal events on the 
low-order Legendre polynomial Kir~ moments. This is 
shown by the shaded histogram in Fig. 1341 which rep- 
resents the difference between the dashed curve in this 
figure and the solid curve in Fig. 120( b). 

A complete representation of the highly localized 
Z(4430)~ distribution in Fig. [3^1 requires the use of Leg- 
endre polynomials to order more than 30, and so our 
low-order representation of the Kn~ angular structure is 



unable to do this. The shaded histogram does reach a 
maximum at about the Z(4430)~ mass value, but corre- 
sponds to a width of ~ 160 MeV, which is almost four 
times larger than the input signal value. 

A fit to the distribution of Fig. [M] with the normal- 
ization of the new Kir~ background, and the normaliza- 
tion, mass and width of the Z(4430) _ signal, free yields 
™z(4430)- = 44 33 ±3 MeV/c 2 , r z(4430) - = 34 ±12 MeV, 
and iVz(44 30 )- = 1402 ± 315 events so that the width is 
reduced to ~ 75%, and the signal to ~ 94%, of the input 
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value, while the mass is essentially unchanged. The solid 
curve in Fig. [M] represents the fit result, and the dot- 
ted curve shows the reduced level of Kir~ background. 
If the Z(4430) _ width is fixed to 45 MeV, we obtain 
iV Z (44 30 )- = 1558 ± 220 events. This is consistent with 
the input value, and presumably more properly reflects 
the effect of the statistical uncertainties in the underlying 
data distribution. We therefore use this to estimate the 
magnitude of the systematic signal reduction factor, and 
so obtain a value of ~ 90%. 





\|r(2S)rc 



(GeV/c 2 ) 



FIG. 34: The m 



i/>(2S)ti 



distribution of Fig. I33f a) for the data 
and MC samples combined as described in the text. The 
dashed curve represents the Kn~ background distribution, 
obtained as described in the text, and the shaded histogram 
represents the difference between this dashed curve and the 
solid curve of Fig. I20f bh i.e. it represents the impact of the 
MC Z(4430)~ signal. The solid curve is the result of a fit 
using the Kn~ background function and a relativistic S-wave 
BW, and the dotted curve shows the resulting renormalized 
dashed curve. 



FIG. 33: (a): The combined m^,(2S)iT- mass distribution for 
the data of Fig.^b) and the MC Z(4430)~ sample generated 
as described in the text; the shaded histogram represents the 
MC sample; (b) the Ktt~ mass distribution corresponding to 
Fig. I33f a) ; the shaded histogram represents the reflection of 
the MC Z(4430)~ events, and the curve shows the result of 
the fit. 

This is a direct demonstration in support of the state- 
ment in the Belle letter [5J that a narrow peak in the 
if)Tt~ mass distribution cannot be generated by only 
P-, and D-wave amplitudes in the Kn~~ system, and in- 
dicates that our use of low-order Legendre polynomials in 
creating our Kix~~ background could lead to an approxi- 
mately 10% systematic reduction of a narrow Z(4430)~ 
signal of the reported magnitude. 

C. Branching fractions 

In Table IVI1 we summarize the branching fraction val- 
ues and their 95% c.l. upper limits, obtained for the 
individual B decay modes studied in the present analysis 
by repeating the fits of Figs. I3"l7a) and l3"TTdh but with 
Z(4430)~ mass and width fixed to the central values ob- 
tained by Belle. The errors quoted are statistical, and 



the upper limits were obtained using these values. For 
the i/j(2S)it~ modes, the branching-fraction and upper- 
limit values have been increased by 10% in order to take 
account of possible reduction of Z(4430)~ signal size, as 
described in Sec. IXIBl 

The branching fraction for the decay mode B° — > 
Z(U30)-K+, Z(4430)- -> ip(2S)n- from Belle is (4.1 ± 

1.0 ± 1.4) x 10~ 5 , to be compared to our upper limit of 

3.1 x 10~ 5 at 95% c.l. 



XII. SUMMARY AND CONCLUSIONS 

We have searched for evidence supporting the existence 
of the i?(4430)~ in the tpTT~ mass distributions result- 
ing from the decays B~'° — > -07r~A' O:+ in a large data 
sample recorded by the BaBAR detector at the PEP-II 
e + e~ collider at SLAC. Since the relevant Dalitz plots are 
dominated by mass and angular distribution structures 
in the Ktt~ system, we decided to investigate the extent 
to which the reflections of these features might describe 
the associated tpir~ mass distributions. To this end, we 
obtained a detailed description of the mass and angu- 
lar structures of the Kir~ system based on the expected 
underlying S-, P-, and D-w&ve Ktt~ amplitude contri- 
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TABLE VI: The Z(4430)~ signal size, the branching fraction value, and its 95% c.l. upper limit for each decay mode; the 
errors quoted are statistical, and the upper limits were obtained using these values; the Z(4430) - mass and width have been 
fixed to the central values obtained by Belle [jj. 



Decay mode 




Z(4430)~ signal 


Branching fraction 
(xlO" 5 ) 


Upper limit 
(xl0~ 5 at 95% c.l.) 


B- -> Z(U30)-R",Z{U30)' 


— ► J/ifin 


-17 ± 140 


-0.1 ±0.8 


1.5 


B" -f ^(4430)-^: + ,Z(4430)- 


— ► J/ip-K~ 


-670 ± 203 


-1.2 ±0.4 


0.4 


B- -» Z{A430)- K°,Z (U30)~ 


-» ij}{2S)ir~ 


148 ± 117 


2.0 ± 1.7 


4.7 


B° -> Z(4430)~iv + ,Z(4430)- 




415 ± 170 


1.9 ±0.8 


3.1 



butions, and in the process even found evidence of D — F 
wave interference in the J/ip decay modes. The frac- 
tional S-, P-, and D-wave intensity contributions to the 
corrected Kn~ mass distributions for B~ — > J/ipir~Kg 
and £>° — > J / K + were found to be the same within 
error, as are the branching fraction values after all cor- 
rections, and so we combined these data in our analysis 
of the J/ipir~ mass distributions. We observed similar 
features for the corresponding ijj(2S) decay modes, even 
though the analysis sample is ~ 6.6 times smaller than 
that for J/ip (Table HJ), and we combined the charged 
and neutral B-meson samples for the ip(2S) analysis also. 

We next investigated the ipir~ mass distributions on 
the basis of our detailed analysis of the Kir~ system. 
We used a MC generator to create large event samples 
for B~'° — ► ijjir~K a ' + with Ktt~ mass distribution gen- 
erated according to the overall fit function obtained from 
the corrected data, but with a uniform cos 9 k distribu- 
tion. The cos 9k dependence was then modulated using 
normalized Legendre polynomial moments whose values 
were obtained from our corrected data by linear interpo- 
lation. 

The total corrected J /i\m~ mass distribution 
(Fig. H0I a )) is well described by this Kn~ back- 
ground, whose form can be seen more clearly in 
Fig. l2"2TaL The residuals (Fig. l2"CTc)) show no evidence 
of a Z(4430)~ signal, and this is true also for the 
various regions of Kir~ mass shown in Fig. [MJa)-(e), 
Fig. HHta),(c), and Fig. |2T[a),(c). When we fit the 
data using a function which allows the presence of a 
Z(4430)~ signal, we obtain only negative Z(4430)~ 
signal intensities (Fig. EUaHc), Table|V|. We find this 
to be the case also for the B~ and B° modes separately, 
and summarize the corresponding branching fraction 
upper limits in Table IVT1 for Z(4430)~ mass and width 
fixed at the central values obtained by Belle [5J . 

We conclude that there is no evidence to support the 
existence of a narrow resonant structure in the J/ipTr~ 
mass distributions for our data on the decay modes 



B°'- -> J/4nr-K +fi . 

The corresponding corrected ip(2S)ir~ distributions 
of Figs. EJflM.td), Figs. MS-® and Figs. ESfl-Cj), 
Figs.l^jfjTb) and (d), and Figs. l2"?T b) and (d) likewise show 
no clear evidence of a narrow signal at the Z(4430)~ mass 
position. 

We have directly compared our uncorrected ip(2S)ir~ 
data with K* veto, to those from Ref. Q; there is no 
evidence of statistically significant difference (Fig. EW c')'). 

In order to quantify our .Z(4430)~ production rate es- 
timates, we fit our total corrected ipnx" mass distribu- 
tions of Figs. l2T)f a) and l2T)f bL using the Kir~ back- 
ground shapes shown in these figures, together with a 
Z(4430)~ line shape. For J/tpir~ we obtain a negative 
signal, while for ip(2S)n~ we obtain a 2.7 standard devia- 
tion signal with fitted width consistent with the value ob- 
tained by Belle, but with central mass value ~ 43 MeV/c 2 
higher than that reported by Belle, which corresponds 
to a ±4.7 standard deviation difference. These fit re- 
sults are shown in Figs. [3lTa) and [3TT dL respectively, 
and are summarized in Table |Vl We repeated these 
fits for the individual decay modes of Eqs.([T])-(rj| with 
Z(4430)~ mass and width fixed to the central values re- 
ported by Belle, and obtained the branching fraction and 
upper limit values summarized in Table IVT1 In particular, 
we find a branching fraction upper limit for the process 
B(B° -> Z(U30)-K+,Z- -► i;(2S)TT-) < 3.1 x 1CT 5 
at 95% c.l., and a corresponding value for the reaction 
B(B- -> Z(M30)~K°,Z- -> ^( 25 >~) < 4 - 7 x 10 ~ 5 
at 95% c.l. We conclude that our analyses provide no 
significant evidence for the existence of the Z(4430)~. 

It will be of great interest to see whether or not the 
Z(4430) _ is confirmed by a future analysis based upon 
a significantly larger data sample than is available at 
present. 
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APPENDIX A: THE DALITZ PLOT 
EFFICIENCY-CORRECTION PROCEDURE 

The efficiency is obtained using samples of Monte Carlo 
events corresponding to the decay processes of Eqs. (1)- 
(4) generated uniformly over the final state Dalitz plot. 

In general, the phase space volume element in the 
Dalitz plot corresponding to the decay B — > ipirK is given 
by: 



dp ~ d(m K7r ) ■ d(m% n ) , 



(Al) 



where m Kn - (m^ n - ) is the invariant mass of the Ktt~ 
(-07T - ) system. 

However, when the efficiency is studied in such recti- 
linear area elements, those elements at the plot boundary 
are partially outside the plot, and this leads to a rather 
cumbersome efficiency treatment. The phase space vol- 
ume element of Eq. (|A1|) may be transformed to 



dp' 



V 



m K Tt d(m Kn ) d(cos9 K ) , (A2) 



i.e. 



dp' ~ p ■ qd(mKit)d(co§Qjz) ■ 



(A3) 



where p is the momentum of the t/j daughter of the B 
in the B rest frame, and q is the momentum of the K 
in the rest frame of the Ktt~ system. This expression is 
such that the phase space density is uniform in cos 6 k at 
a given value of uik-k- ■ 



The range of cos9k is [-1)1], and that of tuk-k- is from 
threshold to mg — m^, so that the resultant "Dalitz Plot" 
is rectangular in shape, with the factor p ■ q representing 
the Jacobian of the variable transformation. A plot of 
this kind can then be used readily to study efficiency be- 
havior over the entire phase space region without the 
problems incurred at the boundary of a conventional 
Dalitz plot (see, for example, Appendix B of Ref. [221]). 

The reconstruction efficiency calculated using the 
Monte Carlo simulated events is parametrized as a func- 
tion of rrijcn- and cos 9k, and then used to correct 
the data by weighting each event by the inverse of its 
parametrized efficiency value. For a given mass inter- 
val / = [m K7r - , rn-K-K- + dmic„-], Iti N be the number 
of generated events, and let n reco , represent the number 
of reconstructed events. The generated cos 9k distribu- 
tion is flat, but in general efficiency effects will cause 
the reconstructed cos 9k distribution to have structure. 
Writing the angular distribution in terms of appropri- 
ately normalized Legendre polynomials, 



dN 



d cos 9 k 



N(P )P (cos9k) 



(A4) 



and, 



dn r 



d cos 9 k 



,^2(Pi)Pi(coa9 K ) (A5) 



where the normalizations are such that, 

J P i (cos6 K )P j (cosO K )d(cos6 K ) = 8 ij , (A6) 

where P t = V2^Y° , and Y° is a spherical harmonic func- 
tion. The value of L is obtained empirically. 

Using this orthogonality condition, the coefficients in 
the expansion are obtained from 

(Pi) =— ^ P J (cos9 K )p^d(cos9 K ), (A7) 

nreco J -1 d COS Ok 

where the integral is given, to a good approximation for 
a large enough MC sample, by Y^h=i" Pj (cos 9k { )- The 
index i runs over the reconstructed events in mass in- 
terval /, such that n reco {Pj} ~ T^i=i° E'j(cos9K i ), and 
the effect of efficiency loss on the angular distribution is 
represented through these coefficients. The absolute effi- 
ciency, calculated as a function of cos 9k and ttik-k-, in 
mass interval /, is then given by 



i?(cos 9k, rriKiT- ) 



With 



N(P )P (cos9 K ) 



En 



N 



(A8) 
(A9) 
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and 



Ei = 2 n reC p(p 3 ) = 2 grr^(cos^) ; (A10) 



N 



N 



for a large enough sample (note that the factor 2 enters 
since (Pq)P (cos 9 K ) = 1/2), Eq. ([A"8]) becomes, 

E(cos9 K ,m K - K -) = E Q 
+Ei P l (cos 9 K ) + ... +E L P L {cos9 K ) . (All) 

The mean value of the Pj(cos0Ki), with i = 1, ...,n rec0 , 
corresponding to mass interval /, is written as (Pj). The 
r.m.s. deviation of the Pj (cos k-) w.r.t. (Pj), a, is given 
by 



(^(cos^,)-^)) 2 

2-^ n — 1 



(A12) 



The error on the mean is then S(Pj) = -^f — and from 

Eq. (Km . 



S(Pj) 



' Etr (p J (cose Kt )-(p 3 )y 

Tlrecoijlreco 1) 



1 



E"reco ( P j( COBe K i )) 2 
i=l n rrr .„ 



- (Pj) 2 



Tl re 



(A13) 



by: 



The uncertainty in the parameter Eq — " r "° is given 



S(Eo) = E a 



N 



(A14) 



and the uncertainty in the coefficient Ej is given by: 



mass intervals from threshold to the maximum value ac- 
cessible. As shown in Sec. I VII Fig. the Ktt~ mass 
dependence of the average efficiency parameter, Eq, de- 
pends on the decay mode of the ip involved, and so 
is obtained by using the MC sample for that particu- 
lar decay mode. The angular dependence represented 
by Ei, Ei,.. .etc. does not depend on the individual ip 
decay mode, and so these coefficients are calculated by 
combining the MC samples for the individual ip modes. 
For the B meson decay processes of Eqs. (p}- (gj), the 
main features of the angular dependence of the efficiency 
are very similar, and so we present the results only for 
B° -> iP(2S)tt- K+ by way of illustration. 

Simulated MC events are subjected to the same re- 
construction and event-selection procedures as those ap- 
plied to the data. For the process of Eq. the cos 9k 
distributions for the surviving MC events are shown for 
each Ktt~ mass interval in Fig. [351 We chose a small 
interval size (0.02) in order to investigate the signifi- 
cant decrease in efficiency observed for cos 9k ~ +1 and 
0.720 < niKrr- < 0.920 GeV/c 2 (Figs. EEI c )-(i)) and 

for cos 9k 1 and 0.970 < m K7r - < 1.270 GeV/c 2 

(Figs. l55T h)- (m) ) . Representation of such localized losses 
requires the use of high-order Legendre polynomials; we 
find that L = 12 yields a satisfactory description, as 
demonstrated by the curves in Fig. [351 The Kir~ mass 
dependence of the resulting values of Ei - En is shown 
in Fig. 1361 and is parametrized in each case by the fifth- 
order polynomial curve shown. These parameterizations, 
together with those describing the Kix~ mass dependence 
of Eq for the individual tfj decay modes (Sec. I VI I) enable 
us to calculate the efficiency at any point in the relevant 
rectangular Dalitz plot, and hence for each event in the 
corresponding data sample. We then assign to each event 
a weight given by the inverse of this efficiency value, and 
by using this weight are able to create efficiency-corrected 
distributions. 



6(Ej 



N 



E£r (^i(cos^)) 2 + (ESC ° ( T os9 " ))r 



(A15) 

For each of the processes represented by Eqs.(l)-(4), 
the efficiency analysis is carried out in 50 MeV/c 2 Ktt~ 



As discussed in Sec. IVI[ the efficiency loss for cos 9k ~ 
+1 is due to the failure to reconstruct low momentum 
charged pions in the laboratory frame, while that for 
cos 9k ~ — 1 is due to the similar loss of low momen- 
tum kaons (Fig. [5]). 



[1] M. Gell-Mann, Baryons And Mesons," Phys. Lett. 8, 214 
(1964). 

[2] G. Trilling, J. Phys. G 33, 1019 (2006). 
[3] C. Amsler ei al. [Particle Data Group], Phys. Lett. B667, 
1 (2008). 

[4] L. Maiani et al, Phys. Rev. D 71, 014028 (2005); L. 
Maiani et al, Phys. Rev. D 72, 031502(R) (2005); L. 
Maiani, A.D. Polosa, and V. Riquer, Phys. Rev. Lett. 



99, 182003 (2007). 

[5] S.-K. Choi et al. [Belle Collaboration], Phys. Rev. Lett. 

100, 142001 (2008). 

[6] The use of charge conjugate reactions is implied through- 
out this paper. 

[7] L. Maiani, A.D. Polosa, and V. Riquer, New Journal of 

Physics, 10, 073004 (2008). 
[8] M. Karliner and H. J. Lipkin, Phys. Lett. B 638, 221 



3G 



30 



20 



10 




100 



O 50 
O 

(73 
•>—> 

c 

>150 



100 



50 



150 



100 



50 



0.620 <m(Kjl)< 0.670 
GeV/c 2 



0.870 < m(Kjl) < 0.920 
GeV/c 2 




1.120 <m(Kjl)< 1.170 
GeV/c 2 




(P) 



-1 -0.5 0.5 1 

cos6 K 



60-0 



40 



20 




100 



50 



150 



100 



50 



150 



100 



50 








).670 < m(K7t) < 0.726 
GeV/c 2 



0.920 < m(K7t) < 0.970 
GeV/c 2 



(g) 



1.170 <m(Kjt)< 1.220 
GeV/c 




i I ■ ■ ■ ■ I ■ ■ ■ i I ■ 



1.420 <m(K7t)< 1.470 
GeV/c 2 




1 -0.5 0.5 1 

cos6 K 



80- 




150 



■ i ■■ 1 1 1 ■■ ■■ i ■ 



0.720 < m(K7l) < 0.77( 
GeV/c 2 



0.970 < m(K7t) < 1.02C 
GeV/c 




L ■ ■ ■ ■ I ■ ■ ■ ■ I ■ ■ ■ ■ I ■ ■ ■ ■ _1 



1.220 < m(Kjl) < 1.27( 
GeV/c 




(m) 



i| "" l ""h 



1.470 <m(K7t)< 1.52C 
GeV/c 2 




-1 -0.5 0.5 1 

cos6 K 



100 



50 




100 



50 



150 



100 



50 



■ i ■■ ■■ i ■■ 1 1 i ■ 



■0.770 < m(K7t) < 0.82f 
GeV/c 2 



1.020 <m(K7t)< 1.07C 
GeVA 




1.270 <m(K7l)< 1.32C 
GeV/c 




(n) 



1.520 <m(K7t)< 1.57C 
GeVA 2 




1 -0.5 0.5 1 

cos6 K 



100 




150 b 



100 



50 



150 



100 



■ i 1 1 ■■ i ■■ ■■ i ■ 



0.820 < m(Kjl) < 0.87C 
I GeV/c 2 



1.070 <m(Kjt)< 1.12C 
GeV/c 



CD 



'1.320 <m(Kjt)< 1.37C 
I GeV/c 2 




(O) 







-1 -0.5 0.5 1 

cosG., 



FIG. 35: The cos 9k distributions in 50 MeV/c 2 Kir mass intervals for the decay mode B° — > ip(2S)ir K + . The points 
represent the data, and the curves show the functions calculated from the moments. 
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FIG. 36: The Kit mass dependence of the coefficients Ei through En for the decay mode B° — » jp(2S)ir K + . The points 
show the calculated values, and the curves result from fits to the data using a fifth-order polynomial. 
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